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ABSTRACT 


Results  of  an  extensive  series  of  repeated  load  triaxial 
tests  on  Dublin  Brown  Glacial  Till  are  reported.  Charts 
are  presented  relating  the  relevant  parameters  of  the 
till  to  dry  density  and  moisture  content.  The  effect  of 
stress  history  is  studied.  The  results  indicate  that 
the  Kirwan  Snaith  nomograph  should  be  modified  and 
suggest  that  the  subgrade  strain  design  criterion  should 
vary  with  the  stiffness  of  the  soil. 

A  recently  developed  low  cost  pneumatic  repeated  loading 
test  apparatus  is  described. 

Skempton's  A  and  B  pore  pressure  parameters  have  been 
used  to  interpret  the  response  of  samples  under  repeated 
loading. 

A  method  of  determining  a  stress  dependent  resilient 
modulus  using  charts  prepared  from  insitu  tests  on  a 
subgrade  is  described. 

Results  of  Laboratory  Pavement  Simulator  tests  on  a 
subgrade  material  are  compared  with  those  obtained  from 
undisturbed  samples  tested  under  repeated  loading  conditions, 
and  also  with  those  obtained  from  a  nomograph  prepared  from 
the  tria.iial  testing  program. 

The  modifications  that  were  made  to  the  computer  program 
DEFPAV  to  include  horizontally  layered  anisotropy  are 
described. 

Initial  results  from  a  testing  programme  to  investigate  the 
properties  affecting  the  behaviour  of  natural  gravels  in 
road  pavements  are  presented. 
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Two  computer  subroutines  for  optimisation  of  road 
pavement  design  are  described. 
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CHAPTER  1 


INTRODUCTION 


This  Final  Technical  Report  for  the  two-year  grant  (Grant  No. 
DAER0-76-6-065)  covers  the  research  carried  out  during  the 
period  September  1976  to  December  1978. 

This  is  a  continuation  of  the  work  carried  out  by  the 
Pavement  Research  Group  at  Trinity  College,  Dublin,  Ireland 
during  the  past  twelve  years.  It  has  previously  been  funded 
by  E.R.O.,  with  five  separate  annual  contracts,  followed  by 
one  three  year  Grant  (No.  DA-ERO- 59 1-73-60017) .  Three  Annual 
Reports  were  provided  for  the  latter  grant,  the  last  being 
dated  April  1976. 

The  object  of  the  recent  research  has  been  to  work  towards  a 
simple  and  rational  method  of  design  for  flexible  pavements 
by  means  of  a  better  understanding  of  the  behaviour  of  the 
constituent  materials,  with  possible  verification  in  the  field. 

An  extensive  series  of  triaxial  tests  was  carried  out  on 
Dublin  glacial  till  in  order  to  study  its  behaviour  under 
repeated  loading,  with  particular  emphasis  on  the  variation 
of  the  relevant  properties  with  dry  density  and  moisture 
content.  The  effect  of  stress  history  was  also  studied.  The 
results  of  this  work  indicate  that  the  nomograph  proposed  in 
the  Report  of  April  1975  requires  certain  modification. 

It  is  also  suggested  that  the  resilient  subgrade  strain 
criterion,  which  is  frequently  used  to  limit  permanent 
deformation  in  the  subgrade,  be  varied  with  the  stiffness 
of  the  soil.  Limited  test  results  indicate  that  the  permanent 
deformation  behaviour  of  a  soil  depends  on  the  degree  of  strain 
of  the  samples. 


1 


In  order  to  carry  out  this  program  of  triaxial  tests  a 
simple  low  cost  pneumatic  repeated  loading  test  apparatus 
has  been  developed,  along  with  a  ram  load  cell  and  lateral 
strain  indicator. 

Skempton's  A  and  B  pore  pressure  parameters  have  been  used 
to  interpret  the  resilient  vertical  and  lateral  response  of 
a  sample  under  repeated  loading. 

A  method  of  determining  the  variation  of  resilient  modulus 
with  stress  level  from  an  insitu  test  on  a  subgrade  using 
prepared  charts  is  given. 

Good  comparison  was  obtained  between  the  stress  dependent 
resilient  modulus  determined  from  the  back  analysis  of 
tests  on  a  glacial  till  subgrade  in  the  Laboratory  Pavement 
Simulator  with  that  recorded  in  triaxial  tests  on  undisturbed 
samples  and  also  with  that  interpolated  from  the  chart 
prepared  from  the  triaxial  testing  programme . 

The  T.C.D.  Finite  Element  Program  has  been  modified  to 
include  horizontally  layered  anisotropy. 

The  initial  results  of  a  series  of  repeated  loading  tests 
on  gravel  pavements  in  the  Laboratory  Pavement  Simulator 
are  reported. 

The  computer  subroutines,  SELPAV  and  OPTPAV,  are  presented 
for  optimisation  in  road  pavement  design.  The  former  is 
incorporated  into  the  proposed  design  system  DEFPAV  and  is 
a  modification  of  that  described  in  our  previous  report. 

The  latter  is  an  independent  package  which  considers  the 
cost  implications  of  a  wide  range  of  design  alternatives, 
including  stage  construction. 

Attempts  made  to  verify  our  design  methods  by  monitoring  the 
performance  of  a  section  of  the  Stillorgan  Dual  Carriageway 
outside  Dublin  were  frustrated  by  vandalism,  construction 
difficulties  and  unsuitable  equipment  for  the  field  use.  It 
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is  considered  impracticable  to  get  successful  field 
verification  of  our  work  with  the  present  resources  of 
the  Pavement  Research  Group.  A  proposal  for  funds  for 
this  work  has  been  sent  to  the  European  Economic  Community 
but  no  reply  has  been  received  to  date. 
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CHAPTER  2 


REPEATED  LOADING  APPARATUS 


2 . 1  Introduction 

The  initial  series  of  repeated  loading  triaxial  tests  were 

carried  out  on  the  mechanically  operated  rig  described  by 

Glynn  (1)  but  with  the  method  of  applying  cell  pressure 

altered  to  the  pneumatic/hydraulic  system  shown  in  Fig.  2.1. 

This  alteration  allowed  a  cyclic  cell  pressure  to  be 

superimposed  on  a  constant  initial  confining  pressure.  The 

system  was  designed  to  allow  a  back  pressure  in  a  sample, 

if  required.  The  mechanical  section  of  the  apparatus  has 

been  in  use  in  Trinity  College  since  1967  and  was  beginning 

to  show  its  age.  Considerable  skill  was  required  by  the 

operator  to  get  the  correct  axial  load  and  the  permanent 

deformation  compensator  required  regular  inspection. 

Furthermore,  the  maximum  axial  stress  which  could  be  applied 

to  a  lOOmm  diameter  sample  was  limited  to  approximately 
2 

120kN/m  .  It  was  therefore  decided  to  develop  and  construct 
a  new  twin  cell  repeated  loading  rig  for  the  series  of  tests 
required  for  our  material  characterisation  program. 

Barksdale  (2)  and  others  have  shown  that  a  passing  wheel 
applied  an  approximately  sinusoidal  or  triangular  type  load 
pulse  to  a  subgrade  element  under  a  road  pavement  of  duration 
of  the  order  of  0.2  -  0.4  secs.  Apart  from  the  rotation  of 
the  axis  of  the  principal  stresses  which  occurs  under  a 
passing  wheel,  the  insitu  stresses  on  a  soil  element  can  be 
modelled  in  a  triaxial  cell  by  an  axial  stress  pulse  and  a 
confining  cell  pressure  pulse  of  similar  waveform  and 
duration  superimposed  on  a  static  confining  pressure  which 
represents  the  insitu  overburden  pressure.  The  load  in  a 
triaxial  cell  must  be  applied  without  impact  and  completely 
removed  between  load  applications  to  allow  the  sample  recover 
for  simulated  rest  periods. 

A  servo-controlled  electro-hydraulic  system  would  give  the 
required  loading  but  was  ruled  out  because  of  its  cost. 
Pneumatic  or  pneumatic/hydraulic  systems  are  considerably 
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cheaper  and  easier  to  construct.  Several  such  apparatii 
which  have  been  reported  in  the  technical  literature,  were 
tried  out  but  were  eventually  discarded  either  because  they 
did  not  meet  our  specific  requirements  or  because  they 
included  equipment  not  readily  available  in  Ireland  or  the 
British  Isles. 

2.2.  New  Repeated  Loading  Testing  Rig 

A  pneumatic  loading  rig  was  therefore  designed  and  built 
at  Trinity  College,  Dublin,  to  meet  our  specific  requirements. 
Details  of  this  rig  are  given  on  Fig.  2.2  the  particular 
features  of  which  are  the  following 

o 

(1)  Air  supply  with  normal  working  pressures  of  700  kN/m  . 

(2)  Axial  load  applied  through  a  70mm  dia  piston  capable  of 

2 

applying  300  kN/m  to  a  102mm  dia  sample  with  an  air 

2 

pressure  of  700  kN/m  . 

(3)  Inlet  air  supply  to  piston  taken  from  an  air  reservoir, 
pressure  of  which  is  adjusted  by  an  air  regulator, 
through  an  on/off  solenoid  which  is  controlled  by  a 
timer  unit,  and  through  a  needle  valve  which  controls 
the  air  supply. 

(4)  Exhaust  air  supply  from  piston  taken  through  an  on/off 
solenoid  which  is  controlled  by  a  timer  and  interval 
delay  relay,  and  through  a  needle  valve. 

(5)  Low  pressure  air  supply  taken  through  a  non  return  valve 
to  the  loading  piston  to  bring  loading  ram  in  contact 
with  the  sample  before  each  load  pulse. 

(6)  Water  reservoir  located  above  test  rig.  This  gives  a 

2 

hydrostatic  cell  pressure  of  14  kN/m  when  the  reservoir 
is  open  to  the  atmosphere. 

(7)  Inlet  air  supply  to  reservoir  taken  through  a  timer 
controlled  solenoid  and  the  exhaust  also  through  a  timer 
controlled  solenoid  and  through  a  pressure  relief  valve 
if  it  is  required  to  increase  the  constant  cell  pressure 
above  that  mentioned  in  (6)  above. 

(8)  Timing  relay  (Tl)  Electromatic  asymmetrical  recycler  with 
time  range  of  0.15-3  secs  on  both  on  and  off  time. 
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(9)  Interval  timer  (T2)  which  is  reset  by  a  relay 
working  from  the  "on  time"  of  the  times  relay  (Tl) . 
Friction  in  the  water  time  delays  the  cell  pressure 
pulse  in  the  triaxial  cell  and  consequently  the  load 
pulse  must  be  delayed  to  get  both  waveforms  in  phase. 

(10)  Delay  on  operation  of  timing  relay  (T3)  worked  from 
(T2)  to  control  the  duration  of  the  axial  load  pulse. 

(11)  Interval  times  (T4)  to  delay  closing  the  exhaust 
solenoid  of  the  loading  piston.  The  circuit  of  this 
unit  which  was  designed  and  constructed  at  T.C.D.  is 
shown  in  Fig.  2.3. 

The  sequence  of  operation  is  as  follows: 

(a)  Timer  relay  Tl  closes  solenoid  A  and  opens  B  to 
pulse  the  air  pressure  in  the  water  reservoir  for 
the  desired  pulse  duration  which  is  set  by  Dial  X. 

The  start  of  this  load  pulse  also  closes  the  relay 
which  in  turn  resets  the  interval  times  T2.  At  the 
end  of  the  desired  period  the  solenoid  A  is  opened 
and  B  closed,  thus  relieving  the  air  pressure. 

(b)  After  the  time  interval  set  on  timer  T2,  the  circuit 
to  the  delay  timer  T3  is  opened  which  in  turn  closes 
solenoid  C  and  opens  solenoid  D,  also  triggering  the 
interval  timer  T4,  which  pulses  the  axial  load  for 
the  length  of  time  set  on  T3. 

(c)  After  the  desired  time  set  on  T4 ,  solenoid  C,  which 
opens  to  relieve  the  axial  load,  is  closed  thus 
allowing  the  low  cell  pressure  supply  to  bring  the 
ram  into  contact  with  the  sample. 

(d)  The  sequence  is  repeated  after  the  period  of  time 
set  on  dial  Y  of  the  timer  unit. 

The  initial  calculation  and  the  development  tests  indicated 

that  the  pneumatic  system  could  be  easily  controlled  to 
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give  a  triangular  or  approximately  sinusoidal  loading  but 
gave  an  unsatisfactory  approximation  of  a  reversed  sine 
waveform  on  the  unloading  phase  which  resulted  in  a 
sawtooth  waveform.  An  automatic  system  was  designed  to 
overcome  this  and  hopefully  give  an  approximately  sinusoidal 
unloading  phase.  This  has  not  been  constructed  to  date  as 
the  response  of  the  system  when  the  sample  is  included 
appears  to  give  a  satisfactory  waveform  on  both  phases. 
Typical  axial  and  cell  pressure  loading  pulses  are  shown 
in  Fig.  2.4 

2.3  Load  and  Deflection  Measurements 


The  methods  for  deflection  and  load  measurements  on  the 
old  mechanically  operated  loading  rig  were  as  described 
by  Glynn  (1) ,  with  the  exception  that  the  axial  load  was 
measured  on  a  load  cell  incorporated  in  the  triaxial 
loading  ram  inside  the  cell.  Initially  it  was  intended 
to  measure  the  load  using  strain  gauges  on  the  inside  of 
a  hollowed  out  portion  of  the  ram.  This  proved  difficult 
to  construct  and  lacked  sensitivity.  The  sensitivity  was 
increased  and  the  construction  simplified  by  using  the 
module  sensor  as  shown  in  Fig.  2.5.  This  load  cell 
generally  gave  satisfactory  performance  but  the  following 
points  were  noted: 

(a)  Sealing  off  the  cell  pressure  from  the  sensor  sometimes 
gave  problems,  particularly  at  high  pressure.  The 
sensor  was  insensitive  to  the  normal  range  cyclic  cell 
pressure. 

(b)  The  increase  in  sensitivity  of  the  load  cell  was  not 
as  great  as  was  expected.  This  may  have  been  caused 
by  elasticity  in  the  system  and  could  probably  be 
eliminated  by  careful  design  and  construction. 

Load  measurements  in  the  new  R.L.T.  rig  were  measured  on 
a  Proving  Ring  outside  the  cell  which  incorporated  an 


7 


LVDT  and  also  on  a  strain  gauged  miniature  proving  ring 
constructed  from  the  loading  ram  inside  the  cell.  The 
latter  type  of  load  cell  was  noted  during  a  visit  to 
Nottingham  University  but  in  order  to  enable  it  to  be 
used  under  water  the  strain  gauges  and  leads  were 
waterproofed.  In  Nottingham,  it  is  used  in  a  non- conductive 
oil  environment. 

The  resilient  and  permanent  deflections  of  the  sample  were 
measured  by  a  twin  LVDT  located  on  top  of  the  cell  as 
shown  in  Fig.  2.6.  Permanent  deflection  was  recorded  by 
the  LVDT  and  checked  by  readings  to  the  top  of  nut  A  by 
micrometers. 

2 . 4  Lateral  Deflection  Measurements 

The  initial  series  of  triaxial  tests  carried  out  in  the 
old  mechanical  rig  did  not  allow  the  lateral  deflection  of 
samples  to  be  measured. 

Problems  were  encountered  in  developing  a  simple  reliable 
system  for  the  measuring  to  small  lateral  deflections  of 
samples  surrounded  by  water.  The  strain  band  as  described 
by  El-Ruwayih (3)  when  constructed  at  T.C.D.  was  not 
sufficiently  sensitive  and  was  affected  by  the  cycled  cell 
pressure.  No  miniature  LVDT  was  readily  available  for  use 
under  water. 

Satisfactory  results  have  recently  been  achieved  with  the 

lateral  strain  indicator  (LSI)  shown  in  Fig.  2.7.  This 

apparatus  is  sensitive,  appears  to  be  unaffected  by  cell 

pressure  and  to  date  has  been  reliable.  The  pressrre  exerted 

2 

by  this  L.S.I.  is  of  the  order  of  0.04kN/m  per  0.01mm 

deflection.  This  means  that  for  lateral  deflection  of  0.2mom 

which  would  correspond  to  vertical  strains  of  the  order  of 

0.004  mm/mm  the  pressure  on  the  side  of  a  sample  at  the 

2 

supports  would  be  approximately  0.8kN/m  .  This  is  not 
therefore  considered  an  important  factor. 
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CHAPTER  3 


REPEATED  LOAD  TRIAXIAL  TESTS  ON  A  GLACIAL  TILL 


3.1  Introduction 


The  subgrade  forms  a  major  element  in  a  road  structure, 
yet  comparatively  little  research  has  been  carried  out  on 
the  behaviour  of  subgrade  materials.  It  is  generally 
agreed  that  methods  for  analysing  road  pavements  have 
reached  a  satisfactory  level  of  development  and  that  further 
work  should  concentrate  on  the  choice  of  suitable  input 
parameters  and  on  verification  by  monitoring  insitu 
performance.  It  is  therefore  of  prime  importance  to 
develop  a  better  understanding  of  the  factors  affecting  the 
behaviour  of  subgrade  materials  under  repeated  loading  and 
under  various  stress  histories.  Hopefully  from  the  knowledge 
obtained,  parameters  necessary  for  input  into  the  design 
systems  can  then  be  related  to  simple  laboratory  or  insitu 
tests. 

Kirwan  and  Snaith  (4)  suggested  that  the  resilient  modulus 
(Mf)  of  a  soil  could  be  related  to  its  Relative  Compaction 
(R  )  and  to  its  Relative  Moisture  Content  (R  ) ,  by  means  of 
a  nomograph  where 


Relative  Compaction  =  - ^situ.dr-X  density - 

Modified  AASHO  maximum  dry  density 


Relative  Moisture  content 


_ Insitu  moisture  content 

Modified  AASHO  optimum  moisture 
content 


The  nomograph  that  was  produced  to  relate  these  parameters 
was  based  on  the  limited  results  that  were  available  to 
them  at  the  time.  It  was  decided  to  initiate  a  program  of 
repeated  load  triaxial  testing  of  soils  at  a  wide  range  of 
dry  densities  and  moisture  contents  in  order  to  investigate 
further  the  Kirwan  Snaith  nomograph,  to  study  the  effect  of 
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stress  history  on  these  parameters  and  to  see  if  the  Rc 
and  concept  could  be  extended  to  the  creep  behaviour. 

3 . 2  Testing  Procedures 

A  glacial  till,  commonly  called  Dublin  Brown  Boulder  Clay, 
was  used  throughout  the  testing  program.  This  forms  the 
predominant  subgrade  in  Ireland.  A  large  quantity  of  the 
soil  was  sieved  through  a  5mm  sieve,  thoroughly  mixed  and 
stored  in  bins  for  testing  at  a  later  date.  The  properties 
of  this  soil  are  given  in  Fig  3.1.  Similar  soil  but 
sometimes  with  different  maximum  stone  size  was  used  by 
Snaith  (5),  Glynn  (1)  and  Casey  (6). 

The  old  twin-cell  repeated  load  triaxial  test  rig  as 
described  by  Glynn  (1)  with  the  modifications  mentioned 
in  Chapter  2  of  this  report  was  used  for  the  series  of 
tests  carried  out  to  date  on  the  material  characterisation 
program.  The  new  rig  is  now  operational  and  the  testing 
program  has  been  moved  over  to  that  equipment. 

The  soil  was  compacted  by  the  method  recommended  by 

Snaith  (5) .  Electro-osmosis  was  sometimes  used  to  bring 

the  samples  close  to  saturation.  The  initial  confining 

stress  (cuj  was  held  constant  throughout  the  tests  at 
2  J 

14  kN/m  with  the  cycled  cell  pressure  (  A03  )  in  phase  with 
the  axial  load  at  h  (  Ao1  -  Ao3  )  unless  specific  tests  were 
made  to  study  the  effect  of  cycled  cell  pressure.  The 
sequence  of  loading  was  varied  to  study  the  effects  of 
various  stress  histories. 

Every  effort  was  made  to  duplicate  each  test.  Nevertheless, 
more  often  than  not  the  inevitable  small  differences  in  the 
dry  density  vd  and  moisture  content  were  sufficient  to 
significantly  affect  the  results.  An  example  of  the  effect 
of  a  small  variation  of  moisture  content  on  the  permanent 
strain  is  shown  in  Fig.  3.2.  Where  samples  were  more  or 
less  the  same,  the  results  from  each  were  consistent,  as 
shown  in  Fig.  3.3.  The  test  results  from  the  different 
cells  were  also  consistent. 
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3 . 3  Modulus  of  Resilience  Versus  Dry  Density  and  Moisture  Content 


The  variation  of  the  modulus  of  resilience  (M  )  with 
and  m/c  for  an  applied  axial  stress  of  20kN/ni  is  shown 
in  Fig.  3.4.  These  results  form  a  consistent  trend  and 
as  can  be  seen,  agree  well  with  those  recorded  by  Snaith  (5) 
which  are  also  plotted  for  comparison. 

One  noticeable  feature  is  the  absence  of  any  local  peak 
value  for  Mr  near  mod  ASSHO  optimum  density  and  moisture 
content  which  would  be  predicted  from  the  Kirwan  and  Snaith 
nomograph.  Contours  of  equal  interpreted  from  the  present 
tests  are  shown.  Replotting  the  results  of  Seed  et  al  (7)  in 
a  similar  form  (Fig.  3.5)  gives  the  same  type  of  contour 
diagram.  It  is  therefore  suggested  that  part  of  the  Kirwan 
Snaith  nomograph  requires  modification.  Contours  of  equal  Mr 
plotted  against  Rc  and  Rw  for  the  brown  glacial  till  are 
shown  in  Fig.  3.6.  The  curves  for  Mod  AASHO  and  Proctor 
Compaction  are  superimposed  on  the  tentative  contour  diagram 
shown  in  Fig.  3.4.  The  maximum  resilient  modulus  was  not  at 
the  optimum  moisture  and  maximum  dry  density  for  these 
compactive  efforts.  A  similar  finding  can  be  interpreted 
from  the  results  of  Seed  et  al  (7).  The  resilient  modulus 
plotted  against  moisture  content  for  given  compactive  efforts 
for  the  T.C.D.  and  Seed  et  al  tests  shown  in  Fig  3.7 
illustrate  the  point.  This  result  is  contrary  to  the  findings 
of  Majedzadeh  and  Guirguis  (8) ,  who  reported  that  for  a  given 
compactive  effort,  the  maximum  resilient  modulus  occurred  at 
optimum  conditions  of  dry  density  and  moisture  content. 

3. 4  Effect  of  Incremental  and  Decremental  Loading  on  the 
Modulus  of  Resilience. 


The  change  in  the  resilient  modulus  with  increases  or 
decreases  in  the  applied  axial  load  and  cell  pressure  are 
illustrated  by  some  typical  results  shown  in  Fig  3.8.  Soft 
samples  showed  little  change  in  Mr  with  increases  in 
deviator  stress,  medium  stiffness  samples  decreased  in 
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modulus  with  increasing  stress  level  while  the  results  of 
the  tests  on  stiff  samples  showed  an  initial  increase  in  M  , 
followed  by  a  decrease.  The  initial  increase  in  Mr  exhibited 
by  the  stiff  soils  was  not  expected.  It  may  have  been  caused 
by  the  seating  of  the  top  cap.  However,  further  tests  are 
being  carried  out  to  investigate  this. 

The  behaviour  of  the  soft  and  medium  strength  stiff  soils  can 
be  explained  by  considering  the  resilient  deformation  as  being 
dependent  on  the  structure  of  the  soil  and  on  the  orientation 
of  the  particles. 

The  particles  of  a  soft  soil  reorientate  themselves  under  low 
stresses  and  any  decrease  in  modulus  with  increase  in  stress 
is  due  solely  to  the  non-linear  elasticity  of  the  soil  system. 
An  increase  in  applied  stress  for  a  medium  strength  stiff  soil 
can  gradually  break  down  the  orientation  of  particles  and 
possibly  lead  to  an  increase  in  compaction.  Therefore  the 
change  in  Mr  reflects  the  non-linear  elasticity  of  the  soil 
structure  as  well  as  the  work  hardening  due  to  compaction. 
Subsequent  unloading  and  loading  at  constant  permanent  strain 
would  be  expected  to  give  lower  moduli  values  should  the 
former  process  predominate,  as  the  structure  would  be  broken 
down  to  some  extent.  The  test  results  shown  in  Fig.  3.9. 

2 

confirm  this  for  the  medium  stiff  soil  (Mr =  50  to  80  MN/m  ) . 
Some  investigators  have  found  an  increase  in  modulus  of 
resilience  for  some  soils  at  high  deviator  stresses.  Under 
these  conditions  the  work  hardening  effect  would  predominate. 

Soil  compacted  at  optimum  modified  AASHO  moisture  content  and 
dry  density  would  be  classified  as  medium  stiffness  soils, 
therefore  the  stress  dependent  modulus  would  be  affected  by 
the  stress  history.  It  is  suggested  that  the  modulus-stress 
relationship  for  the  loading,  unloading  curve  be  used  for 
design  on  the  basis  that  it  is  the  weaker  condition  and  also 
the  condition  which  would  exist  after  an  initial  settling-in 
period  of  a  subgrade.  Expressing  the  stress  deoending 
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dependent  resilient  modulus  in  the  form: 


M 


r 


the  test  results  for  the  unloading,  reloading  curves  gave 
an  approximate  value  of  of  0.14. 


The  sine  function  relationship  developed  previously  by 
Kirwan  and  Snaith  (4)  did  not  fit  these  curves  but  most  of 
the  present  tests  were  incrementally  loaded  as  opposed  to 
individual  samples  being  tested  under  different  loads. 


3.5  Permanent  Strain 


In  order  to  study  the  variation  of  permanent  strain  (e  ) 

P 

with  dry  density  yd  and  moisture  content  (m/c)  it  was 
necessary  to  model  the  behaviour  of  each  sample  by  a  creep 
function  to  compare  parameters.  The  permanent  strain  for 
9  tests  are  plotted  in  Fig.  3.10  against  deviator  stress 
for  40,000  cycles.  Owing  to  the  lack  of  adequate  information 
on  the  creep  behaviour  of  soils,  permanent  deformation  of 
subgrade  materials  is  generally  controlled  in  road  design  by 
ensuring  that  the  maximum  resilient  subgrade  strain  does  not 
exceed  values  which  have  been  determined  from  the  back  analyses 
of  existing  road  structures  which  have  behaved  satisfactorily. 
The  relationship  suggested  by  Brown  et  al  (9)  for  British 
soils  gives  a  limiting  subgrade  strain  of  1.0  x  10  ‘  for  a 
design  life  of  40,000  cycles.  Working  back  from  the  measured 
values  of  the  resilient  modulus  from  the  test  results  the 
resulting  permanent  strains  are  shown  in  Fig.  3.10.  Thus 
the  working  range  is  less  than  1.0%  permanent  axial  strain. 
Below  this  value  the  permanent  strain  results  can  be 
approximated  to  a  straight  line  ^  =  fn(N,  m/c)*?^. 
where  N  =  number  of  cycles.  This  function  fn(N,^,  m/c) 
corresponds  to  the  creep  compliance  used  by  Glynn  (1)  the 
value  of  which  (at  a  given  number  of  cycles)  can  be  used 
to  compare  the  permanent  deformation  of  soils  of  different 
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dry  densities  and  moisture  contents.  Permanent  strain  is 

generally  expressed  in  the  form  =  fn  (N)  a  for 

P 

constant  values  of  moisture  content  and  dry  density. 

Glynn  (1)  on  his  tests  on  Dublin  brown  boulder  clay  (Fig. 

3.11)  used  a  value  of  b  =  1  thus  (assuming  linear  variation 
with  stress)  but  incorporated  the  concept  of  a  threshold 
stress.  Applying  the  condition  that  the  applied  load  must 
result  in  strains  below  those  allowed  by  the  subgrade  strain 
criterion,  Glynn's  test  results  can  be  approximated  by  a 
straight  line  if  the  threshold  stress  is  ignored.  Kirwan  and 
Snaith  (10)  argued  that  a  value  of  b  =  1.75  was  more  correct 
but  their  results  are  more  difficult  to  interpret  as  separate 
samples  were  used  at  each  stress  level  with  corresponding 
small  differences  in  yd  and  m/c.  Nevertheless  the 
assumption  of  linear  behaviour  appears  to  be  reasonable  within 
the  normal  working  range  of  stress  and  deformation. 

3.6  Creep  Compliance  vs  Dry  Density  and  Moisture  Content 

The  creep  compliance  at  40,000  load  repetitions  for  the  tests 
to  date  are  shown  in  Fig  3.12,  plotted  against  dry  density  and 
moisture  content.  The  results  form  a  consistent  trend  and 
illustrate  the  dramatic  increase  in  permanent  deformation  that 
can  occur  with  an  increase  in  moisture  content. 

3. 7  Subgrade  Strain  Criterion 

The  subgrade  resilient  strain  criterion  would  be  expected  to 
give  a  relatively  constant  permanent  deformation  if,  as  it  is 
generally  assumed,  it  controls  the  permanent  deformation  in 
the  subgrade.  Most  of  these  empirically  derived  relationships 
depend  solely  on  the  number  of  load  repetitions  and  therefore, 
like  that  of  Brown  et  al  which  is  plotted  on  Fig.  3.13,  would 

give  much  less  permanent  deformation  for  the  stiff  samples 

2 

(0.2%  for  M  =  75MN/m")  than  for  the  softer  samples  (1%  for 
Mr  =  20MN/in  )  .  Barker  et  al  (11)  from  a  back  analysis  of 
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airfield  runways,  suggested  that  the  limiting  design  subgrade 
strain  be  varied  with  the  stiffness  of  the  soil.  Their 
relationship  would  allow  an  increase  in  the  resilient  strain 
for  stiffer  soils,  which  is  in  agreement  with  the  test  results. 
The  tests  to  date  therefore  suggest  that  the  subgrade  strain 
criterion  should  be  varied  with  the  stiffness  of  the  soil. 

3. 8  Permanent  Strain  and  Stress  History 

The  limited  results  to  date  suggest  that  the  permanent  strain 

is  strain  dependant.  This  is  shown  in  Fig.  3.14  where  two 

identical  samples  were  incrementally  loaded  with  different 

initial  stress  levels.  It  can  be  seen  that  once  the  samples 

reached  the  same  permanent  strain  their  behaviour  was  almost 

2 

identical.  A  sample  unloaded  from  a  stress  level  of  80kN/m 
and  subsequently  unloaded  and  reloaded  at  various  lower  stress 
levels  suffered  no  measurable  increase  in  permanent  deformation 
as  can  be  seen  in  Fig.  3.15. 

The  test  results  indicate  that  when  a  certain  strain  is 
reached  as  a  result  of  a  given  number  of  repetitions  of  load, 
at  a  particular  stress,  then  additional  repetitions  of  load 
at  a  lower  stress  will  not  result  in  any  further  increase  in 
strain  provided  that  the  strain  attained  at  the  higher  stress 
level  exceeds  the  maximum  possible  strain  at  the  lower  stress 
level.  This  behaviour  agrees  with  a  "time-hardening  model" 
for  permanent  strain  suggested  by  Monismith  et  al  (12)  but 
here  the  term  strain  dependent  is  considered  more  suitable. 

This  point  is  illustrated  in  Fig.  3.16.  The  cumulative 
permanent  deformation  of  a  spectrum  of  loads  cannot  therefore 
be  added  together  to  give  the  total  deformation,  as  the 
sequence  of  loading  would  affect  the  result.  This  is  a  point 
that  is  being  investigated  further  at  T.C.D. 

3.9  Permanent  Strain  vs  Resilient  Strains 


Barker  (11A)  found  a  relationship  between  the  permanent  and 
resilient  strains  from  published  results  of  repeated  load 

i 
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tests  from  various  investigators.  The  results  of  the  present 
series  of  tests  plotted  in  Pig.  3.17  to  both  a  natural  and 
log-log  scale  tend  to  confirm  this  finding.  The  results  of 
Snaith  (5)  fall  within  the  same  zone. 
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Fig.  3.5 


RELATIVE  MOISTURE  CONTENT 
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Fig,  3.10 


w 


PERMANENT  STRAIN  Vs  DEVIATOR  STRESS  AT  40,000  CYCLES 


Fig.  3.11 
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CHAPTER  4 


Study  of  effect  of  degree  of  saturation  on  Resilient 
Modulus  and  Poisson's  Ratio. 

The  lateral  deflection  of  samples  under  repeated  loading 
in  the  triaxial  apparatus  was  not  measured  in  the  initial 
series  of  tests  because  of  the  difficulty  in  developing 
measuring  equipment  of  sufficient  sensitivity.  In  order 
to  simulate  the  permanent  deformation  of  the  soil  insitu, 
the  axial  loads  and  lateral  pressure  were  cycled  in  phase. 

The  measured  vertical  deformation  was  therefore  a  combination 
of  deformation  under  axial  load  and  volume  compressibility 
under  the  cell  pressure  e.g. 

e  =  (Aox  -  Act3)  Aq3  where  e  =  Vertical  Strain 

v  —  M  *  2  K  ^ 

r  r  Mr  =  Resilient  Modulus 

or  Kr  =  Bulk  Modulus 

v  =  Poissons  Ratio 

e  =  (Act i  -  2vAa3) 
v  M 

r 

assuming  that  the  material  behaved  elastically. 

Poisson's  Ratio  can  be  estimated  by  closing  off  the  cycled 
cell  pressure  for  one  cycle,  measuring  the  Resilient  Modulus, 
and  then  substituting  back  into  the  equations  above.  Given 
that  Poisson's  ratio  varies  between  0.1  and  0.5  normally, 
it  is  clear  that  the  measured  Mr  and  therefore  the  Mr  com¬ 
puted  from  tests  in  which  the  cell  pressure  is  cycled,  can  be 
in  error  by 

Measured  Modulus  =  A°i  “  Aa3  =  Agi  ~  Ag3 

e  Aa i  -  2vAa3 

V  w 


if  v  =  0.5  i.e.  the  material  is  saturated,  the  Resilient 


n 


Modulus  computed  from  the  vertical  deflection  with  the  cycled 
cell  pressure  would  be  similar  to  that  computed  without  the 
cycled  cell  pressure.  The  error  between  the  measured  value 
of  Mr  and  the  computed  value  of  is  of  the  order  of  25%  for 

v  =  0.2  and  a  cycled  cell  pressure  of  H(AOj  -  A o3) 

It  is  of  interest  to  study  the  expected  variation  with  and 
m/c.  To  do  this,  use  is  made  of  the  Skempton's  A  and  B  pore 
pressure  parameters  (13) .  The  material  is  assumed  to  be 
elastic  (a  staged  approached  where  the  permanent  deformation 
if  negligible) ,  linear,  isotropic  and  that  the  degree  of 
saturation  is  sufficiently  high  such  that  surface  tension  can 
be  ignored. 


With  reference  to  Bishop  and  Henkel  (14),  the  change  in  pore 
pressure  due  to  applied  total  stresses  Aa3  and  Aa3  is 

Au  =  B{ Aa 3  +  A  ( Aa  3  -  Aa3)} 
for  isotropic  elastic  behaviour  A  =  1/3 
.’.  Au  =  B{Ao3  +  j(Ao3  -  Acr3)} 

The  vertical  strain  is  therefore: 


e  = 
v 


AO; 


-  Ao3 

T1 


(l-2v  1 ) 

E1 


{AO  3  ~  B  ( Aa  3  + 


(Acti  -  Ap3) 


where  E'  =  Elastic  modulus  of  soil  skeleton  in  effective 
stresses 

v1  =  Poisson's  ratio  of  soil  skeleton  in  effective 
stresses . 


Taking  Aa3  =  Y (Aa i  -  Aa3) 
e  =  (Acti  -  Aa 3 )  {  1  +  (1-2vM  (Y  -  B  (Y  +  i)  )} 

^  fl  J 
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The  resilient  modulus  measured  in  a  triaxial  sample  with  no 
cycled  cell  pressure  is  therefore: 


{i-|(l-2v1)  }  {  l-|(l-2v1)  } 


with  the  cell  pressure  cycles  at  *5  (A Oj  -  Ao3)  as  it  was  for 
the  initial  tests.  The  deviator  stress  divided  by  vertical 
strain  is 

Aoj  -  Aa3  E' 

ev  {  1  +  (l-2v 1  )  (0.5-0.833B)  } 

The  variation  of  the  ratio,  Aoj  -  Ao3  is  shown  in  Fig.  4.1. 


assuming  that  v'  =  0.2  as  got  from  Calladine  (15)  for  a  soil 
with  Plasticity  Index  of  14. 

This  again  shows  that  the  error  in  the  estimation  of  Resilient 
Modulus  is  negligible  at  100%  saturation  but  can  be  expected 
to  increase  at  lesser  saturation  levels.  A  knowledge  of  the 
variation  of  the  B  pore  pressure  parameter  with  Yd  and  m/c  would 
therefore  give  an  indication  of  the  degree  of  error. 

Extending  this  analysis  to  study  the  measured  resilient  Poisson's 
Ratio  under  repeated  loading,  the  lateral  strain  is 

er  =  ^  (A011  +  Ao3l) 

E'  E' 

=  jjj,  {Ac3  -  Au  -  vl(Aaj  +  Ac3  -  2Au)  } 
neglecting  the  cycles  cell  pressure 

er  =  {-Au-v 1 (Ao 1  -  2Au)  } 
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but 


AU  =  B  AOi 
3 


er  -  {f  +  v1  (1-|£)  } 

Hi  J  J 


V  = 

r 


|®L  {|  +  v1  (1~|^)  } 


U-B  ( l-2v  1 )  } 


|  +  v1  (1--^) 


1  -  |(l-2vJ) 


.  . . (2) 


It  can  be  shown  that  equation  (2)  is  also  true  if  the  cell 
pressure  is  cycled. 

With  B  =  1  i.e.  saturated  vr  =0.5 

as  expected,  regardless  of  the  Poisson's  Ratio  of  the  soil 
structure. 

With  3  =  0.5  v  =  0.166+0.67 

r  _ 

0.834  +  0.332 

Thus,  although  the  Poisson's  Ratio  of  the  structure  could  be 
0.2,  the  resilient  Poisson's  Ratio  of  the  partially  saturated 
material  could  be  much  higher. 

With  B  =  0  vr  =  v1  obviously 

With  =  0.25,  a  low  value  B  =  0.5  vr  =  0.364 


The  variation  in  the  resilient  Poisson's  ratio  against  the 
B  pore  pressure  parameter  is  illustrated  in  Fig.  4.2.  Thus 
the  measured  Poisson's  Ratio  in  a  repeated  load  triaxial 
apparatus  would  be  expected  to  be  higher  than  that  of  the 
soil  skeleton  owing  to  the  build-up  of  positive  pore  pressure 
during  loading. 


Equation  (1)  is  independent  of  the  cycled  cell  pressure,  thus 
the  Poisson's  Ratio  measured  with  a  constant  confining  pressure 
should  be  similar  to  that  measured  with  a  cycled  confining 
pressure.  Allen  and  Thompson  (15a)  recorded  noticable 
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different  values  for  Poisson's  Ratio  for  these  two  types  of 
tests  which  suggests  that  either  the  soil  structure  changed, 
or  that  there  were  inaccuracies  in  the  measuring  equipment, 
or  possibly  anisotropic  behaviour  of  the  samples. 

Preliminary  triaxial  tests  have  confirmed  that  the  resilient 
vertical  deflection  is  unaffected  by  cell  pressure  for 
saturated  sampled  but  that,  for  unsaturated  samples,  the  cell 
pressure  had  a  slightly  greater  effect  than  predicted  by  pore 
pressure  parameter  analysis.  Previous  investigations  have 
found  th^t  the  effective  stress  elastic  modulus  of  some  soils 
is  dependent  on  the  confining  pressure.  Work  is  being  carried 
out  to  see  if  this  is  significant,  particularly  as  high  soil 
suction  or  negative  pore  water  pressures  are  involved  in 
partially  saturated  soils.  Considering  the  modulus  of 
resilience  in  terms  of  a  constant  effective  stress  modulus 
(Equation  1)  implies  that  a  reduction  of  the  B  pore  pressure 
parameter,  i.e.  drying  out  of  the  soil,  reduces  the  resilient 
modulus.  This  is  not  true  in  practice.  Thus,  while  this 
form  of  analysis  can  usefully  study  the  behaviour  of  a  single 
sample  under  varying  stress  conditions,  as  yet  it  cannot  be 
extended  to  predict  the  effect  of  changing  moisture  content 
or  dry  density.  Information  on  soil  suction  variation  may 
bridge  this  gap.  This  is  being  investigated. 
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CHAPTER  5 


METHOD  OF  DETERMINING  RESILIENT  MODULI  FROM  DEFLECTION 
MEASUREMENTS 


5 . 1  Introduction 

For  pavement  design  procedures  based  on  a  computer  analysis 
of  stresses,  strains  and  deflections  in  a  multi-layered  system, 
an  estimate  of  the  elastic  moduli  of  each  constituent  material 
is  required  as  input  to  program.  The  most  difficult  material 
to  characterise  in  this  way  is  the  subgrade,  since  its  modulus 
depends  on  moisture  content,  dry  density,  stress  level, 
structure  and  stress  history. 

Subgrade  material  can  be  tested  in  repeated  load  triaxial 
testing  machines  so  as  to  determine  the  relative  effects  of 
the  above  mentioned  variables  on  the  dynamic  response  of  the 
soil.  While  this  test  is  necessary  to  determine  the  suscept¬ 
ibility  of  a  soil  to  creep  deformation  and  to  provide  an 
estimate  of  the  variation  in  modulus  with  traffic  loading  and 
changes  in  environmental  conditions,  it  has  its  limitations 
in  yielding  a  suitable  design  modulus  for  a  subgrade,  since 
it  is  impossible  to  model  the  true  in-situ  stresses  and  in- 
situ  soil  structure  in  a  triaxial  type  machine. 

In  recent  years,  the  results  from  many  full-scale  instrumented 
test  tracks  have  been  used  to  provide  verification  of  the 
predictive  capabilities  of  the  various  theoretical  analysis 
procedures.  In  general,  it  has  been  found  that  the  greatest 
difficulty  has  arisen  in  selecting  suitable  characteristics 
for  the  subgrade  material. 

A  simple  method  has  been  devised  to  interpret  the  results  of 
an  in-situ  deflection  test  on  the  subgrade  to  obtain  an 
estimate  of  its  stress  dependent  modulus.  It  is  based  on  the 
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principle  that  the  shape  of  the  deflection  bowl  is  a 
function  of  the  rate  of  variation  of  modulus  with  stress 
right  through  the  soil  mass. 

Other  investigators  (Leger  (16),  Grant  and  Walker  (17),  among 
others)  have  used  a  radius  of  curvature  as  a  measure  of 
deflection  bowl  profile,  to  obtain  information  on  base  and 
subgrade  moduli  in  roads  with  thin  surfacings.  Also,  deflection 
ratios  have  been  used  to  obtain  estimates  of  moduli  in  a  multi¬ 
layered  elastic  system.  Valkering  (18) ,  by  assuming  a  subgrade 
with  a  constant  modulus,  overlain  by  a  granular  base  course 
having  a  constant  modulus  which  is  a  fixed  ratio  of  that  of 
the  subgrade  (Dormon  and  Metcalf  (19)  )  has  produced  a  chart 
relating  the  modulus  of  the  surface  layer  to  the  deflection 
ratio.  This  principle  has  been  extended  for  use  in  the  Shell 
(Claessen  and  Ditmarsch  (20)  )  overlay  design  method. 

Ullidtz  (21)  has  also  used  deflection  ratios  in  his  overlay 
design  method  and  allowed  the  non-linear  effect  of  the  sub¬ 
grade  to  be  taken  into  account  in  their  interpretation.  Others 
have  used  measured  deflections  at  off-centre  points  to  produce 
a  spreadability  index  which  can  be  related  to  material  moduli 
values . 

Here  the  use  of  deflection  ratios  is  extended  to  allow  the 
stress  dependent  modulus  of  a  subgrade  to  be  determined  with 
the  aid  of  a  non-linear  elastic  computer  analysis.  From  this, 
nomographs  have  been  produced  which  allow  the  stress  dependent 
modulus  to  be  estimated  without  using  the  program.  Also,  the 
method  of  obtaining  the  in-situ  modulus  values  of  the  consituent 
layers  in  a  multi-layered  system  is  explained. 


.  2  Basis  of  Nomograph 

The  principle  behind  the  nomograph  can  best  be  understood  by 
considering  the  finite  element  modelling  of  a  single  layer 
isotropic  material  with  boundary  conditions  and  loading  area 
constant. 
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where  { P }  =  load  vector 

[k]  =  stiffness  matrix 

{ U }  =  deflection  vector 

The  stress  dependent  modulus  of  resilience  can  be  modelled  by 
a  function  of  the  stress  in  an  element  e.g.  Mr  =  Cfn(o(j) . 

where  C  =  constant  a,  =  deviator  stress 

d 

The  individual  terms  in  the  stiffness  matrix  are  multiplied  by 
the  modulus  of  resilience  appropriate  to  the  stress  condition 
in  the  element.  For  a  single  layer  the  constant  in  equation 
(2)  above  can  be  taken  out  of  the  stiffness  matrix  such  that 
[k]  =  C  [k-] 

where  [K']  is  the  adjusted  stiffness  matrix. 

Equation  (1)  therefore  becomes 

(P }  =  C  [k']  {U  } 

{u}  =  1  [K^  fp}  (3) 

C 

Since  each  term  on  the  right-hand  side  of  equation  (3)  is 
multiplied  by  1/C,  the  relative  magnitude  of  the  terms  in  {U } 
remains  independent  of  C.  The  stress  dependent  modulus  of 
resilience  of  soils  is  generally  represented  by  an  exponential 
function  e.g. 

Mr  =  Cod  "Y  where  C  and  y  are  soil  constants . (4) 

a.  is  the  deviator  stress 
d 

The  deviator  stress  is  proportional  to  applied  load  for  a 
linear  elastic  system.  Maher  (22)  has  shown,  using  the  finite 
element  program  DEFPAV,  that  for  a  constant  value  of  y  a 
similar  relationship  exists.  Thus 

i  P, 

-  OC  — ■ 

Go  Po 

For  an  increase  in  applied  load  from  PQ  to  P^,  each  modulus  of 
resilience  in  the  stiffness  matrix  is  increased  by  (P^/Po)  "*Y 

from  equation  (4)  above.  This  has  no  effect  on  the  ratio  of 


deflections  but  decreases  each  term  by  1/(P^/PQ). 

Therefore  (a)  the  relative  magnitude  of  terms  in  the 

deflection  vector  {U}  is  determined  solely 
by  the  value  of  y 

(b)  the  magnitude  of  the  deflection  is  proportional 
t°  p 
cIpT^ 


The  nomograph  is  drawn  up  for  c  =  100  and  P  =  100, 
Given  two  soils  with  modulus  relationships. 

M__ ,  =  C,  a ,  ^  and  M  _  =  C_a  Y 


rl 


'lud 


r2 


'2d 


under  loads  {P^  }  and  {P2  }  respectively,  the  ratio  of  the 
central  deflections  are  from  equation  (5) : 

Pi 

Ul 


_  CiP 


-Y 


u2 


(£«) 

C2  P  2 


(E1) 


Cl  —  C 2  U  2  p  2 


C  1  ^P_l  ^  “  Y 


C2P2 


-Y 


Ui  fPi ,-Y 


The  nomograph  shown  in  Fig.  5.1  shows  the  relationship  between 
Y  and  the  deflection  for  a  range  of  radial  distances,  computed 
from  DEFPAV. 


Also  plotted  on  the  nomograph  is  the  central  deflection  (l^) 
computed  by  DEFPAV,  for  C2  =  100  and  P2  =  100,  for  various 
values  of  y.  The  value  of  for  the  soil  measured  can  there¬ 
fore  be  calculated  from 


Ci 


u  2 

100  x  — 


m 


x 


P 

a 

100 


-Y 


where  u  is  the  measured  central  deflection  and  P 
in  si 

applied  pressure. 


is 


the 


5 . 3  Boundary  Effects  and  Influence  of  Poisson's  Ratio 

Since  the  interpretation  method  is  based  on  measured 
deflections,  it  is  essential  that  the  computer  model,  from  which 
the  nomograph  is  prepared,  is  an  exact  representation  of  the 
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1 


measured  system. 

Particularly  in  the  case  of  the  simulator  tank  results, 
boundary  discontinuities  must  be  placed  in  the  approximate 
positions  in  the  computer  model.  Fig. 5. 2  shows  the  effect  of 
promixity  of  the  vertical  boundary  (nodes  free  to  move 
vertically  but  fixed  horizontially)  on  the  values  of  Dr. 

Beyond  1.5  metres  the  effect  is  negligible. 

Since  the  elastic  behaviour  of  a  subgrade  is  dependent  on  the 
value  of  Poisson's  Ratio  as  well  as  on  that  of  the  Resilient 
Modulus,  it  is  important  to  determine  the  effect  of  this 
parameter  on  the  values  of  the  deflection  ratios. 

Fig. 5. 3  shows  the  variation  in  deflection  ratios  for  Poisson's 
Ratio  varying  from  0.22  to  0.47.  It  can  be  seen  not  to  have 
a  marked  effect,  thus  making  the  use  of  a  constant  value  an 
acceptable  simplification. 

To  illustrate  the  use  of  the  interpretation  methos,  suppose  a 

2 

clay  subgrade  is  loaded  with  a  stress  of  50kN/m  over  a  circular 
area  of  radius  100mm  and  a  deflection  of  0.27mm  is  measured  at 
the  centre  and  0.08mm  at  120mm  from  the  centre.  Then 

°120  =  =  0.296 


Then,  from  Fig. 5.1,  for  0^20  =  0.296,  y 
«=td  =  2*7mm'  Maher  (22)  suggests. 
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o  s  td 


x  p 
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0.85  and 


2.7  x  50 

0.27  x  ( 5 0 )~ s ’ 8 5 
100 


=  227.8 

M  =  277.8  ad"0,85  MN/m2 

•  IT 


51 


5.4  APPLICATION  OF  THE  METHOD  IN  THE  INTERPRETATION  pfr  SOME 
PREVIOUSLY  PUBLISHED  RESEARCH 

Tests  were  performed  at  the  University  of  Nottingham  to 
measure  the  stresses  and  strains  occurring  in  a  typical 
cohesive  subgrade  when  subjected  to  dynamic  loading.  The 
tests  were  performed  in  a  test  pit,  1.52m  deep  by  2.44m 
square  with  the  load  being  applied  to  the  surface  over  a 
circular  area  of  radius  152mm  using  a  pneumatic  loading 
system. 

A  series  of  tests  were  carried  out  (Sparrow  &  Tory (24)  )  at 

2 

varying  contact  pressures  (50-100kN/m  )  and  loading  times. 

Surface  deflections  were  plotted  in  terms  of  normalised  values 

so  that,  for  the  purposes  of  this  investigation,  they  have 

been  converted  to  actual  deflections  corresponding  to  an 

2 

applied  stress  of  69kN/m  (Fig. 5. 4). 

Assuming  a  semi-infinite  elastic  isotropic  soil  with  a  Poisson's 

2 

Ratio  of  0.3  then  a  modulus  of  12.4  MN/m  is  required  to  produce 
a  similar  centreline  deflection  to  the  measured  value.  However, 
at  off-centre  points  it  can  be  seen  that  the  deflection  bowl 
departs  substantially  from  that  measured. 

By  producing  nomographs  in  the  same  form  as  that  in  Fig. 5.1 
but  with  radius  of  152mm  and  boundaries  similar  to  those  in 
the  tests,  then  a  stress  dependent  modulus  can  be  determined 
which  results  in  deflections  which  more  closely  follow  the 
measured  values.  Exact  coincidence  is  not  obtained  due  to  the 
fact  that,  as  reported  by  Sparrow  and  Tory,  the  resilient 
modulus,  calculated  from  in- situ  measurements,  increased 
markedly  with  depth  and  also  the  exact  influence  of  the 
boundaries  is  unknown. 

Approximate  linear  and  exponential  Mp/a^  relationships  for  the 
test  material  are  shown  in  Fig.  5.5. 

A  further  more  thorough  series  of  tests  were  performed  by 
Brown  and  Pell  (25)  in  the  same  test  pit.  The  effect  of  applied 
stress,  loading  area  and  loading  frequency  on  the  measured 
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response  of  the  soil  was  investigated.  It  was  found  that, 
to  satisfactorily  predict  the  subgrade  strains  and  maximum 
surface  deflections,  using  linear  elastic  theory,  a  variety 
of  modulus  values  was  required.  It  was  found  that  strains 
and  deflections  were  also  very  dependant  on  local  values  of 
the  Resilient  Modulus. 

Fig  5.6  shows  the  measured  deflection  points  for  four  of  the 
tests.  A  summary  of  the  test  conditions  is  given  in  Table  5.1. 


Test 

No. 

Normal 

Contact 

Pressure 

(kN/M2) 

Time 
for  load 
pulse 
(sec. ) 

Estimated 

Y 

CP 

S/A 

117 

2.0 

0.38 

67.64 

S/B 

117 

0.1 

0.44 

93.15 

S/C 

52 

0.1 

0.435 

78.10 

S/D 

52 

2.0 

0.33 

56.35 

Table  5.1  Summary  of  Nottingham  Tests 


It  can  be  seen  from  Fig.  5.6  that  fitting  the  Boussinesq 
equation  to  the  centre-line  deflection  causes  a  large  divergence 
from  the  measured  values  with  distance  from  the  loading  area. (25) 

Assuming  an  exponential  Mr/a<3  relationship  and  obtaining 
average  values  for  Y  and  Cp  (Table  5.1)  which  most  closely  fit 
the  deflection  ratio  at  each  of  the  measured  points,  Mr/  od 
curves  as  shown  in  Fig.  5.7  are  obtained.  It  is  seen  that 
these  fall  very  closely  together  and  give  a  more  accurate  repre¬ 
sentation  of  the  measured  deflection  values  (Fig.  5.6)  at  off¬ 
load  points.  It  is  also  seen  that,  when  a  constant  value  of 
modulus  is  chosen,  the  value  required  to  yield  the  correct 
centreline  deflection  in  each  case  is  very  dependent  on  the 
magnitude  of  the  test  conditions. 
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5.5  INTERPRETATION  METHOD  FOR  MULTI -LAYERED  SYSTEMS 


For  a  two-layer  linear  elastic  system  with  fixed  dimensions 
and  subjected  to  a  constant  stress  over  a  fixed  circular  area, 
the  resilient  moduli  of  the  material  can  be  determined  from 
the  shape  and  magnitude  of  the  deflection  bowl.  The  shape  of 
the  deflected  profile  is  a  function  of  the  ratio  of  the  two 
moduli  and  the  maximum  centre-line  deflection  depends  on  the 
absolute  value  of  these  moduli.  Thus  an  interpretation 
nomograph  can  be  prepared  giving  the  relationship  between  the 
deflection  ratio,  the  ratio  of  the  moduli  and  the  magnitude 
of  the  centre-line  deflection,  for  predetermined  values  of 
layer  thickness  and  the  load  radius.  The  applied  stress  can 
be  accounted  for  by  plotting  the  deflection  in  terms  of 
deflection  per  unit  stress  or  more  conveniently  as  deflection 
per  100  kN/M^  of  applied  stress. 

In  the  case  of  a  two- layer  system  in  which  the  lower  layer  is 
a  clay  subgrade  with  a  marked  non-linear  modulus,  the  deflection 
bowl  shape  as  measured  at  the  surface  will  be  substantially 
modified,  particularly  where  the  surface  layer  is  thin  and  of  a 
comparable  modulus  to  that  of  the  subgrade.  To  allow  for  this, 
the  value  of  one  of  the  parameters  in  the  modulus/stress 
relationship  of  the  subgrade  must  be  assumed.  The  nomograph 
can  then  be  prepared  in  terms  of  this  known  parameter  and  the 
unknown  one  estimated  from  the  nomograph  using  the  measured 
deflections . 

Examples  of  both  these  nomographs  and  their  use  are  given  in 
Chapter  8 . 

The  value  of  the  vertical  stress  at  either  of  the  layer  inter¬ 
faces,  if  known,  also  provides  information  on  the  modular 
ratios  and,  in  conjunction  with  the  deflection  ratio,  deflection 
magnitudes  and  layer  thicknesses  would  permit  determination  of 
the  layer  moduli. 
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5 . 6  SUMMARY 


A  simple  method  is  derived  and  explained  for  back-analysing 
measured  deflection  profiles  from  dynamic  tests  on  a  subgrade. 
This  provides  the  magnitude  and  non-linearity  of  the 
resilient  modulus  as  a  function  of  deviator  stress. 

The  technique  is  demonstrated  for  an  assumed  exponential 
modulus  stress  relationship.  However,  the  method  is  not 
restricted  to  this.  Nomographs  can  be  constructed  in  terms 
of  any  desired  function  and  for  any  given  loading  conditions. 

Examples  are  given  which  demonstrate  the  use  of  the  analysis 
procedure  in  providing  consistent  modulus  values  for  a  series 
of  test  results. 
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CHAPTER  6 


Comparison  of  Resilient  Modulus  determined  from  Pavement 
Simulaton  tests  to  those  from  triaxial  tests  and  from  the 
nomograph. 

6.1  INTRODUCTION 


In  order  to  take  the  study  of  the  behaviour  of  glacial  till 
under  repeated  loading  a  stage  beyond  triaxial  testing  and 
nearer  the  insitu  conditions,  two  series  of  tests  were  carried 
out  on  glacial  till  subgrades  in  the  Laboratory  Pavement 
Simulator.  Undisturbed  samples  of  these  subgrades  were  also 
recovered  and  tested  in  the  repeated  loading  triaxial  rig  to 
study  the  effects  of  sampling  procedures  on  the  relevant 
parameters . 

The  Laboratory  Pavement  Simulator  (L.P.S.)  has  previously  been 
described  by  Kirwan  et  al  (26) .  The  resilient  and  permanent 
deflections  were  measured  at  the  centreline  and  at  distances 
of  120mm  and  210mm  from  the  centre  on  an  extended  diameter. 


6 . 2  TESTING  METHODS 


The  tests  were  carried  out  on  the  glacial  till  subgrade  after 
removal  of  the  gravel  sub-base  used  in  the  Gravel  Road  experim¬ 
ents  of  the  Foras  Forbartha  which  are  described  in  Chapter  8 
of  this  report.  The  surface  was  covered  with  wax  immediately 
after  removal  of  the  gravel.  The  glacial  till  was  taken  from 
a  different  area  of  Dublin  (Ballybrack)  than  that  used  in  the 
triaxial  experiments  and  consequently  exhibited  slightly 
different  gradings  and  index  properties  as  can  be  seen  from 
Fig.  6.1. 

The  following  is  a  summary  of  the  testing  program  for  the  two 
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series  of  tests  carried  out  in  the  L.P.S.  to  date.  In  each 

series,  three  of  four  areas  of  the  subgrade  surface  were 

loaded  vertically  under  applied  uniform  pressure  which  varied 
2  2 

from  100kN/m  to  300kN/m  over  a  circular  area  of  204mm 
diameter.  The  resilient  and  permanent  deformations  were 
measured  at  the  centre  and  at  distances  of  120mm  and  210mm 
from  the  centre. 


L.P.S.  Test  Series  A 


Average  m/c  11.1% 

2 

Average  dry  density  19 . 5mg/m 


Insitu  CBR 


25  to  29% 


Test  No. 

Applied 

No.  of  stress 

pressure 

kN/m2 

application 

C8 

70-320 

C9 

73 

1,600 

154 

1,900 

CIO 

238 

10,000 

303 

10,000 

181 

9,000 

Cll 

70-310 

C12 

130 

2,000 

205 

4,000 

C13 

40-330 
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L.P.S.  Test  Series  B 


Average  m/c 
Average  dry  density 
CBR 


11.4% 

1. 98Mg/m' 
15-18% 


Test  No. 

Applied 

No.  of  stress 

pressure 

«v 

application 

kN/ni 

1.  a) 

100 

6,000 

b) 

200 

67,000 

2.  a) 

100 

32,000 

b) 

200 

44,000 

c) 

300 

29,000 

d) 

100-300 

Varied 

3.  a) 

100 

30,000 

b) 

200 

21,450 

c) 

100-200 

Varied 

d) 

300 

28.000 

e) 

100-300 

Varied 

4 ,  a) 

200 

26,000 

b) 

300 

78,000 

c) 

100-300 

Varied 

Undisturbed  sampled  were  taken  following 

each  test,  both 

adjacent  to  and  from 

under  the  test  areas 

.  These  samples 

were  obtained  either  by  jacking  U4  thin  walled  sample  tubes 
(lOOnun  dia,  450mm  in  length)  into  the  soil  or  by  cutting 
blocks  of  the  soil. 


6.3  RESILIENT  RESPONSE 

SERIES  A 

The  results  of  this  test  series  indicated  much  higher  moduli 
than  anticipated.  The  recorded  readings  for  Series  A  must 
therefore  be  questioned.  No  reason  can  be  found  for  such  an 
error  other  than  that  the  tests  were  carried  out  under  super¬ 
vision  as  part  of  a  final  year  project  by  students  who  were 
not  completely  familiar  with  the  equipment.  Further  tests 
are  in  progress  which  hopefully  will  confirm  that  the  Series  A 
results  were  in  error.  The  report  on  this  test  series  given 
below  should  be  read  with  this  in  mind. 


Centreline  deflections  only  were  measured  in  tests  C8  and  C9 . 
The  material  was  assumed  linear  and  the  modulus  of  resilience 
determined  by  Boussinesq  analysis.  This  gave  values  of: 


C8  Mr  =  340MN/m2  at  100kN/m2  dropping  to  283KN/m2  at  300kN/m2 

C9  Mr  =  628MN/m2  at  100kN/m2  dropping  to  501MN/m2  at  300kN/m2 

The  resilient  moduli  for  tests  CIO  and  C13  were  determined 

from  the  deflection  bowl  measurements  using  the  nomograph 

given  in  Fig.  6.2.  This  assumes  a  stress  dependent  modulus 

of  the  form:  =  c  -Y 

r  « 

where  C  and  y  are  material  parameters  and  the  deviator 
stress.  The  results  are  shown  in  Fig.  6.3.  Good  agreement 
was  generally  obtained  for  individual  tests  carried  out  at 
different  loading  pressures  but  it  was  clear  that  the  sub¬ 
grades  for  the  latter  two  tests  C12  and  C13  were  considerably 
softer  than  the  previous  tests. 


The  resilient  modulus  of  the  above  test,  which  comprised  the 

Series  A  tests,  were  surprisingly  high  for  brown  glacial  till 

at  the  measured  moisture  content  and  dry  density.  This  would 

also  explain  the  agreement  between  the  empirical  relationship, 

2 

M  =  10  x  CBR  MN/m  and  the  reported  insitu  values. 
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Repeated  load  triaxial  tests  on  undisturbed  samples  gave  a 

2 

Resilient  Modulus  of  38  to  54MN/m  which  was  almost  independent 

of  the  deviator  stress.  The  lower  modulus  values  were  on 

samples  taken  from  previously  loaded  areas.  These  results  are 

considerably  lower  than  those  recorded  on  the  subgrade  in  the 

L.P.S.  The  constand  confining  pressure  in  the  triaxial  tests 

2 

which  was  normally  14kN/m  to  represent  typical  over-burden 

2 

pressure,  was  increased  to  70kN/m  with  no  marked  effect  on  the 
measured  resilient  response.  The  samples  were  remoulded  and 
recompacted  to  their  original  moisture  content  and  dry  density 
but  again  gave  similar  results.  The  relative  densities  were: 


Sample 

Number 


1 

2 


Dry 

Density 

rr 

Mg/m 


Moisture  Relative 

Content  Density 

Percentage 


Relative 

Moisture 

Content 


1.968  11.1%  0.95  1.11 

1.918  11.0%  0.93  1.10 


The  resilient  modulus  determined  from  the  nomograph  given  in 

Fig.  3.6,  Chapter  3  using  the  above  Yd  and  m/c  9ave  values  of 

2  2 

approximately  60kN/n\  for  a  deviator  stress  of  20kN/m  .  This 

compares  favourably  with  the  values  measured  in  the  triaxial 

tests  but  not  to  those  from  insitu  tests  in  the  L.P.S.  The 

back  analysis  of  the  tests  on  the  Gravel  Pavements  in  the  L.P.S. 

(Chapter  8)  showed  a  subgrade  modulus  of  the  order  of 

80-110MN/m  .  Assuming  that  a  stress  dependent  modulus  for 

the  subgrade  gave  lower  values  as  shown  in  Fig.  6.4.  This 

agrees  favourably  with  the  triaxial  results  and  the  nomograph 

but  not  with  those  obtained  from  tests  on  the  surface  of  the 

subgrade  or  with  the  empirical  relationship  relating  the 

resilient  modulus  to  the  CBR.  It  would  appear  therefore  that 

there  were  recording  errors  in  the  results  of  those  tests  taken 

on  the  subgrade  surface. 

SERIES  B 

The  resilient  moduli  recorded  from  the  second  series  of  tests 


* 
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on  a  subgrade.  Series  B,  were  considerably  lower  than  those 
reported  for  Series  A. 


The  stress  dependent  moduli  (assuming  M  =  C  ad"Y)  for  tests 

3  and  4  are  shown  in  Fig.  6.5.  Test  2,  an  area  previously 

loaded  in  the  gravel  experiments,  gave  an  almost  stress 

2 

independent  modulus  of  between  163MN/m  with  an  applied 

2  2 

pressure  of  lOOkN/m  ,  to  94MN/m  at  an  applied  pressure  of 
300kN/m2 . 

Tests  carried  out  at  different  applied  pressures  but  on  the 
same  test  area  of  the  subgrade  surface  gave  consistent  results. 

A  total  of  twelve  undisturbed  samples  were  taken  from  the 
subgrade  after  the  test,  three  from  under  the  loaded  area,  six 
from  areas  not  previously  loaded  and  three  taken  horizontally 
to  study  anisotropy. 

Two  samples  have  been  tested  to  date,  all  in  the  recently 
developed  Repeated  Loading  Testing  Rig.  The  results  of 
these  two  tests,  shown  in  Fig.  6.6,  show  similar  type  curves 
for  Mr  Vs  as  that  given  by  back  analysis  of  the  insitu  tests. 

Furthermore,  the  results  of  the  two  tests  are  compatible.  The 
triaxial  results  are  slightly  lower  than  the  insitu  values 
which  is  to  be  expected  from  sampling  disturbance.  The  dry 
densities,  moisture  content,  relative  density  and  relative 


moisture 

content 

of  the  samples 

were  as 

follows : 

Sample 

Dry 

Density 

Mg/m3 

Moisture 

Content 

Percentage 

Relative 

Density 

Relative 

Moisture 

Content 

6 

1.974 

11.5% 

0.954 

1.15 

A 

1.989 

11.6% 

0.96 

1.16 

The  resilient  moduli  interpreted  from  the  nomograph  (Fig. 3. 6, 

2  2 

Chapter  3)  are  55MN/m  and  50MN/m  for  the  two  samples.  These 


compare  favourably  with  the  measured  insitu  response.  Further¬ 
more,  back  analysis  of  the  Gravel  tests  which  were  carried  out 

on  this  subgrade  gave  resilient  moduli,  approximately 
2 

45-64  MN/m  ,  which  again  compares  favourably. 


6.4  PERMANENT  DEFORMATION 


The  results  of  the  Series  A  tests  were  discarded  because  of 
possible  recording  errors  which  became  apparent  from  the 
analysis  of  the  resilient  response. 

The  permanent  deformation  against  number  of  load  repetitions 

for  tests  1,  2  and  3  are  shown  in  Fig.  6.7  and  for  test  4,  in 

Fig.  6.8.  The  latter  tests  are  shown  separately  as  the  applied 

2 

load  reduced  from  200kN/m  to  zero  when  run  overnight  (from 
6,000  to  19,000  repetitions).  The  general  shape  of  the  curves 
is  similar  to  those  obtained  from  the  triaxial  tests  reported 
in  Chapter  3. 

The  triaxial  tests  reported  in  Chapter  3  indicated  that  , 

within  the  normal  working  subgrade  stress  range,  the  permanent 

deformation  could  be  taken  as  an  approximately  linear  function 

of  deviator  stress  for  a  given  number  of  repetitions.  This 

implies  that  there  is  a  straight  line  relationship  between  the 

permanent  deformation  in  the  L.P.S.  tests  and  the  applied 

pressure.  These  parameters  are  plotted  on  Fig.  6.9.  Test  2 

did  give  a  straight  line  plot  but  this  area  had  been  loaded 

previously  in  Foras  Forbartha  gravel  experiments.  The  tests 

gave  an  almost  linear  response  for  applied  pressure  of  up  to 

200kN/m3  but  diverged  from  a  straight  line  when  the  applied 

2 

pressure  was  increased  to  300kN/m  .  The  maximum  resilient 

2 

strain  in  the  subgrade  under  300kN/m  pressure  was  of  the  order 
of  3.5  x  10"3.  The  curve  produced  by  Brown  et  al  (9)  implies 
that  the  subgrade  under  this  resilient  strain  should  fail  with 
512  repetitions.  This  load  is  therefore  outside  the  normal 
level  which  would  be  applied  to  the  subgrade. 
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The  permanent  deformation  vs  applied  pressure  can  be 

approximated  to  a  straight  line  for  the  loads  less  than  • 

2 

200kN/m  as  shown  m  Fig.  6.9.  The  slope  of  the  permanent 

deformation  -  deviator  stress  line  obtained  from  the  * 

triaxial  tests  was  called  the  modulus  of  compliance  (Mp)  and 

was  used  to  compare  the  performance  of  different  samples.  The 

line  shown  on  Fig. 6. 9  implies  a  modulus  of  compliance  of 
-5  2 

2  x  10  /kN/m  .  This  compares  favourably  with  that  obtained 
from  the  triaxial  tests  for  soils  of  similar  stiffness. 

The  Relative  Compaction  vs  Relative  Moisture  Content  plot 

2 

shown  on  Fig. 3.6,  Chapter  3  indicated  an  Mp  »  0.00009  kN/m 
for  the  subgrade  tested.  This  is  greater  than  that  recorded 
by  a  factor  of  4.5.  The  contours  of  equal  Mp  change  rapidly 
for  a  small  change  in  moisture  content  in  the  area  of  the 
diagram  in  which  the  insitu  soil  tests  are  plotted.  Never¬ 
theless,  the  reasons  for  this  descrepancy  are  being  investigated. 

f 

6.5  SUMMARY 

Two  series  of  repeated  load  tests  were  carried  out  on 

glacial  till  from  Ballybrack,  Co.  Dublin  in  the  Laboratory 

Pavement  Simulator.  Analysis  of  the  resilient  response  of 

these  tests  indicated  possible  recording  errors  in  the  first 

series  of  tests.  The  stress  dependent  resilient  modulus 

determined  by  the  deflection  bowl  method  described  in  Chapter  5 

on  the  second  series  of  tests  compared  favourably  with  those 

recorded  on  triaxial  tests  on  undisturbed  samples  and  with 

values  interpolated  from  the  chart  given  on  Fig. 3. 6,  Chapter  3. 

2 

The  empirical  relationship  Mr  =  10CBR  MN/m  gave  unsatisfactory 
results . 

While  the  general  trends  of  the  permanent  deformation  complied 

with  the  behaviour  recorded  for  the  triaxial  tests,  the  actual 

deflection  was  approximately  4.5  times  less  than  predicted 

from  the  R  ,  R  chart  given  on  Fig. 3. 13,  Chapter  3. 
c  w 


69 


DEVIATOR  STRESS  0^  KN/M: 


TEST  2 


CHAPTER 


ANISOTROPIC  VERSION  OF  DEFPAV 


Several  investigators  have  reported  different  moduli  of 
resilience  for  samples  taken  vertically  to  those  taken 
horizontally  from  road  pavement  layers  (Shackel  (27)  ) .  This 
is  not  surprising  considering  that  the  method  of  construction 
involves  compaction  in  layers  using  a  roller.  It  is  therefore 
important  that  methods  of  analysis  are  available  to  study  the 
effect  of  the  layered  anisotropy  on  the  road  response  and 
performance.  Such  methods  may  also  allow  more  satisfactory 
modelling  of  fabric  layers  than  those  which  use  transition 
layers  with  isotropic  properties.  The  Trinity  College,  Dublin 
finite  element  program  DEFPAV  has  been  modified  to  include 
horizontally  layered  anisotropic  elastic  material  properties. 

The  formulation  of  the  finite  element  program  DEFPAV  for 
isotropic  material  properties  has  been  described  by  Kirwan  and 
Glynn  (28) .  Modification  to  layered  anisotropic  behaviour 
involves  altering  the  [d]  matrix  which  relates  stresses  in 
individual  elements  to  strains 
{a}  =  [d]  { e } 

This  modification  of  stress-strain  relationship  effects  the 
element  stiffness  matrix 

[k]  =  £  [Bj  [d]  [b]  dv 

where  [bD  relates  strains  to  nodal  displacements  in  local  co¬ 
ordinates. 

The  stresses  are  computed  from  the  relationship 

{a}  =  [d]  [b]  [aT1  {u} 

where  A  relates  nodal  displacement  in  global  co-ordinates 
to  nodal  displacement  in  local  co-ordinates. 


The  derivation  of  the  Cd]  matrix  for  layered  anisotropic 
elastic  materials  assumes  that  the  vertical  direction  will 
be  the  axis  of  symmetry  for  any  anisotropy  which  may  exist. 
This  implies  that  the  governing  parameters  are  independent 
of  the  horizontal  direction.  Some  aspects  of  the  behaviour 
of  anisotropic  soils  have  been  discussed  by  Pickering  (1970) 


DEFPAV  used  a  rectangular  annular  element  as  shown  in  Fig. 7.1. 
The  strains  in  a  section  of  the  element  shown  in  Fig.  7.1  are 
given  by: 


°z  _ 

eT 

0 

vzh  _  “ 

Eh 

vzh  _1 
Eh 

er  c 

"Miz 

Ez 

♦  - 

% 

v  a6 

vhh _ 

% 

e0  = 

-^hz!i 

Ez 

- 

Eh 

+  °9 

*h 

Y 

r  z  « 

1_  Trz 

G  2 

Where 

oz,ez  represent  stress  and  strain  in  the  vertical  direction. 

oz  , ez  represent  stress  and  strain  in  the  radial  direction. 

00 , £0  represent  stress  and  strain  in  the  circumferential 

direction . 

Trz  *  Yrz  represent  stress  and  strain  in  shear. 

Vzh  represents  the  degree  of  vertical  strain  due  to 
horizontal  strain. 

\hz  represents  the  degree  of  horizontal  strain  due  to 
vertical  strain. 


vhh  represents  the  degree  of  horizontal  strain  due  to 
horizontal  strain  in  the  orthogonal  direction. 

G*  represents  the  shear  modulus  in  the  vertical  plane. 


The  input  requirements  of  the  program  are  similar  to  those  of 
the  isotropic  version  however  they  must  be  extended  to  include 
the  terms  illustrated  in  APPENDIX  1. 
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From  energy  considerations 


h  2 

Let  —  =  =r-  =  m  and  ordering  the  matrix  in  the  same 

b  hi 

z  z 

way  as  used  by  Kirwan  and  Glynn  (28),  a  modified  [d]  matrix 
which  relates  stresses  to  strains  in  a  horizontally  layered 
anisotropic  material  is  then  obtained.  This  [d]  matrix  is 
illustrated  in  Fig.  7.2. 

The  [b]  matrix  remains  unaltered  thus  the  stiffness  matrix 
in  local  coordinates  is  given  by: 

[k]  =  ///  [B j"  [d]  [b]  dz .  dr.  r  d0 
Fig.  7.3  shows  a  simplified  form  of  this  stiffness  matrix 
where  the  following  substitutions  have  been  made; 

Ml  bn(l-nvhz)  .  in  (r^/rd 

M2  bn(l+v,  J  (r  .-r.  ) 

tin  j  l 

M3  bn(l+v  )  (r2. -r2.  ) 

hn  j  1 

M4  (bV2 ) n (I-nv  ) Hn (r . /r . ) 

hz  j  l 

M5  fcfn(l+v,)  (r  .-r.  ) 

nh  j  i 

M6  (b3/3)n(l-n\£>  Intr^/r.  ) +bm{  (l+vhh)  ( l-v^nv^  ( r^  -r2. )  /2} 
M7  (b2/2)n(l+vhh)  ( r ^ -r. ) 

M8  (b5'2)n(l+vhh)  (r'.-r2.) 

M9  (bJ'SJnd+v^)  (rj-r.)+bm{  (1+v^)  (l-\^h~2nv^  (r3.-!3.  )/3} 

M10 :  -  (b3/3)  n  ( l+vhh)  ( r2.  -r2. ) +bm{ ( l+vhh)  ( l-v^nv^  ( r\  -r“. )  /4 } 

MU:"  Inb(1+Vbh)  (1_Vhh'2nVhV  (rj-ri)/2 

M12:-  mb(l+v  )  (1-v  -2nv2)  (r3. -r3)/3 
hh  hh  hz  j  i 

M13: ”  mb<1+vhh)  (1-Vhh"2nvhz)  (rj-ri)/2 

M14:-  bn  ( 1+v,  J  (r  . -r .  ) 
hn  j  l 

Ml5 : *  bnV1+Vhh)  (rj-r2i) 

M16:-  (b^lnv^l^Xr.-r.) 

Ml? 5 “  ^V1+Vhh>(rrr2i)/2 

M18 : -  b { 1-v2 )  (r2.-r2  )/2 

nh  j  i 


M19  :  -  bnv  ( l+vuJ  ( r2.  -r2.  )  /2 
hz  tin  j  1 

M20:  -  2bnv,  (1+uJ  (r3.-r3.  )/3 
hz  nn  j  1 

M21 : -  (b2/2)nx|5z(l+\^h)  (r2.-r2. ) /2+m(b2/2)  {  d+vhh)  (l-v^nv^  (r2. -r^)  /2 } 

M22 : -  t^ivh2(l+vhh)  (r3.-r3.)/3+m(b2/2)  {  (1+v^)  (l-vhh-2nv£)  (r3.-r3.)/3 

M2 3  :  -  m(b^2)  (1-v  -2nv2)  (r2.-r2.)/2 
hh  nz  j  l 

M24  :  -  b  ( 1-v 2 )  (r3. -r3.  )/3 

hh  j  l 

M25  :  -  bd-v^l  ( r ^ -r“.  ) /4+m (b^3 )  {  <  1+v^)  ( l-v^^nv^  /2 

The  individual  terms  of  the  matrix  products  [d]  x  [b]  for  the 
anisotropic  material  are  included  in  the  subroutine  FIN6EL  in 
DEFPAV  and  are  illustrated  in  Fig.  7.4.  This  matrix  is  some¬ 
what  different  to  that  used  in  the  isotropic  version  of  DEFPAV, 
as  a  result  modifications  have  been  made  to  the  subroutines 
DATAIN,  FIN8EL,  FIN2EL,  FIN6EL  and  FINSTR  in  the  program 
as  well  as  other  minor  alterations  elsewhere. 


It  is  clear  from  the  [k]  matrix  given  in  Fig.  7.3  that  the 

stiffness  becomes  indeterminate  with 

1-v, .  -2nv,  =  0 

hh  hz 

Pickering  (29)  has  carried  out  an  interesting  study  on  the 
limiting  values  of  elastic  parameters  of  horizontally 
layered  anisotropic  materials,  in  order  to  keep  strain 
energy  positive  the  following  conditions  must  hold: 


a)  Ez  and  E^  must  be  greater  than  zero. 

b)  E  n_  .  _  2 

{  ■=-  hh  ^Vzh  }  {  1  +  v,,  }  must  be  positive. 

£*  hn 

Z 

This  is  similar  to 


1  - 


v,  -  2nv2 
hh  hz 


>  0 


c)  Gj  must  be  positive 

d)  2  (  1  +  v,  )  must  be  positive  i.e.  v  ,>-  1 

hn  hh 


Condition  (b)  can  restrict  what  seem  to  be  reasonable 
values  for  the  elastic  parameters. 


IN  CYLINDRICAL  COORDINATES 


SIMPLIFIED  ANISOTROPIC  ELEMENT  STIFFNESS  MATRIX  FOR  ANNULAR  RECTANGULAR  ELEMENT 
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CHAPTER  8 


GRAVEL  PAVEMENT  TESTS 


8.1  INTRODUCTION 

In  Ireland,  it  has  been  found  that  certain  pit  run  glacial 
gravels  behave  satisfactorily  as  pavement  base  course 
materials  while  others  give  poor  results.  For  gravels 
without  excessive  fines  passing  the  75  B  S  Sieve,  there  is 
often  no  obvious  reason  for  poor  performance.  This  is  of 
particular  importance  in  this  country  where  surface  dressed 
gravel  often  forms  an  initial  part  of  the  stage  construction 
of  trunk  roads.  (The  gravel  base  course  with  only  a  thin 
surface  dressing  is  open  to  traffic  for  periods  of  up  to 
several  years  before  application  of  the  bituminous  layers). 

It  has  been  found  that  grading  tests  on  gravel  do  not  provide 
an  indication  of  future  performance.  Compaction  tests  and 
C.3.R.  tests  are  similarly  inadequate.  An  Foras  Forbartha 
therefore  decided  to  carry  out  a  series  of  repeated  load 
tests  in  our  laboratory  pavement  simulator.  A  base  course  of 
pit  run  Kilpeddar  gravel  placed  on  a  subgrade  of  brown  glacial 
till  from  Ballybrack,  Dublin  has  been  tested  to  date.  The 
laboratory  construction  was  to  simulate,  as  closely  as  possible, 
the  conditions  at  the  road  improvement  scheme  at  Kilpeddar, 
County  Wicklow  which  is  being  monitored  in  the  field.  The 
preliminary  analysis  of  this  work  is  reported  in  "The  use  of 
gravels  in  road  pavements  in  Ireland"  by  S.Davitt,  Research 
Officer,  An  Foras  Forbartha.  The  tests  are  continuing  with 
gravels  from  other  sources. 

8.2  PAVEMENT  SIMULATOR  TESTS 

While  the  problem  of  the  behaviour  of  gravels  is  still  unsolved 
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at  present,  some  interesting  trends  have  been  observed.  A 
brief  outline  of  the  laboratory  tests  is  as  follows: 

Subgrade:  Brown  glacial  till:  1  metre  thick 

Series  A:  Water  Content  11.1%  Dry  Density:  1.95  Mg/m3 

C.B.R.  :  25-29%  thin  bituminous  coating  on  surface  of  subgrade 

Base  Course  Gravel:  Water  Content  5%  -  5.5% 

Thickness  250mm  and  150mm 

Insitu  Density:  2.22Mg/m3 

Grading,  compaction  characteristics  and  C.B.R.  of  gravel  given 
on  Fig.  8.1. 

Moisture  content  of  gravel  after  test  =  3% 

2 

Applied  repeated  load:  552  kN/m 
Radius  of  load  pulse:  102mm 
Duration  of  load  pulse:  0.4  sec. 

Load  frequency:  2  cycles  /  sec. 

Stresses  (Vertical) at  gravel/subgrade  interface  were  measured 
using  load  cells. 

Tests  were  normally  carried  out  up  to  100,000  repetitions  of  load 
application.  Surface  deflections  were  measured  on  centre  line 
and  at  120mm  and  210mm  each  side  of  the  centre  on  a  diametral  line. 

Tests  were  carried  out: 

1)  On  the  gravel  at  placement  water  content 

2)  After  increasing  the  water  content  of  the  gravel  by 
soaking  for  various  periods  of  time. 

3)  With  ponding  of  water  round  the  point  of  application 
of  load,  to  achieve  saturation  of  the  gravel. 

The  following  trends  were  observed  on  the  behaviour  of  the 
whole  pavement: 

1)  'Dry':  After  10,000  loading  cycles,  the  permanent  deformation 
of  the  pavement  showed  a  very  small  increase  with  further 
applications  of  load.  The  total  permanent  deformation  for 
150mm  thickness  of  gravel  was  less  them  3mm  after  100,000 
cycles. 
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2)  After  soaking:  an  increase  was  observed  in  the  resilient 
deflection,  and  also  in  the  initial  permanent  deformation. 
Again,  after  10,000  cycles,  when  no  further  moisture 
entered  the  system,  very  little  further  permanent  defor¬ 
mation  occurred. 

3)  Ponding  of  load  application  area:  The  permanent  defor¬ 
mation  as  expected,  increased  rapidly  due  to  the  loading 
pad  punching  into  the  gravel.  Resilient  modulus  also 
decreased  as  expected. 

It  is  concluded  that  the  unsaturated  gravel  base  had  an 
extremely  high  load  carrying  capacity  under  the  vertical  load 
applied  in  the  Laboratory  Pavement  Simulator.  The  extent  to 
which  the  confining  effect  of  the  walls  influenced  the  results 
has  not  been  ascertained. 

As  expected,  the  increase  in  water  content  towards  saturation 
had  an  adverse  effect  which  emphasised  the  importance  of  both 
the  free  draining  requirements  of  the  gravel  and  the  need  to 
slope  the  formation  to  ensure  rapid  drainage  between  base  and 
subgrade.  It  may  be  advantageous  to  blind  the  surface  of  the 
basecourse  to  limit  ingress  of  water  in  the  event  of  surface 
dressing  breakdown. 

The  subgrade  has  been  tested  insitu  and  in  triaxial  repeated 
load  tests  as  part  of  our  materials  characterisation,  these 
tests  are  reported  and  discussed  in  Chapter  6. 


8.3  INTERPRETATION  OF  RESILIENT  RESPONSE  BY  BACK  ANALYSIS 


A  modulus  interpretation  system  has  been  developed  for  the 
two  layer  pavements  used  in  the  L.P.S.  This  method  is  based 
on  estimation  of  modulus  of  resilience  from  deflection  basin 
measurements  assuming  isotropic,  homogeneous  materials  as 
described  in  Chapter  5. 
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The  stress  dependence  of  the  resilient  modulus  is  therefore 
ignored,  which  is  certainly  incorrect  for  this  gravel. 

Despite  these  assumptions,  this  method  is  helpful  in 
analysing  the  behaviour  of  the  materials.  This  analysis 
has  enabled  the  behaviour  of  the  gravel  to  be  separated  from 
that  of  the  subgrade  and  to  compare  the  calculated  stresses 
with  those  recorded  on  the  embedded  stress  cells. 

Interpretation  nomographs  corresponding  to  the  two  pavement 
thicknesses  are  shown  in  Figs.  8.2  and  8.3.  A  curve  relating 
the  stress  on  the  subgrade  to  the  deflection  ratio  is  included 
for  comparison  with  the  measured  values. 

The  results  of  this  analysis  are  summarised  on  Tables  8.1  to 
8.6.  The  main  trends  which  can  be  observed  from  these  are: 

1)  Dry  (ref  GC1  to  GC4)  :  The  resilient  modulus  of  the 
gravel  layer  increased  considerably  during  the  tests. 

That  of  the  subgrade  appears  to  increase  slightly 
although  this  is  partly  due  to  a  reduction  in  the 
subgrade  stress  level. 

2)  Soaked  (ref  GC5  to  GC6)  :  The  resilient  modulus  of  the 
soaked  gravel  decreased  significantly  with  increase  in  the 
number  of  cycles  while  that  of  the  subgrade  increased 
slightly,  probably  as  a  result  of  compaction  of  the  sub¬ 
grade.  The  increase  in  resilient  deflection  therefore 
took  place  in  the  gravel. 

3)  Ponded  :  The  resilient  modulus  of  the  gravel  beneath  a 
free  water  surface  decreased  by  half  during  the  test, 
presumably  due  to  pore  pressure  build-up.  The  increase 
in  modulus  of  the  subgrade  is  unexplained. 

4)  Reduction  in  thickness  of  gravel  from  250  to  150mm: 

The  modulus  of  the  gravel  increased  during  the  first  test 
and  reached  values  slightly  lower  than  those  obtained 
from  the  thicker  pavement.  This  is  in  keeping  with  the 
trend  noted  by  Klomp  and  Dorman  (30). 
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5)  The  average  ratio,  between  the  modulus  of  the  gravel 
layer  and  that  of  the  subgrade,  was  six  for  the  thick 
layers  and  three  for  the  thinner  layers.  These  are  higher 
ratios  than  those  found  by  the  Shell  investigators. 

61  Using  the  constant  modulus  for  both  subgrade  and  gravel, 
the  calculated  stresses  at  the  top  of  the  subgrade  do 
not  always  compare  favourably  with  those  measured.  Table  8.4. 
But  taking  a  stress  dependent  modulus  for  the  subgrade  of 
the  form  M  =  the  comparison  between  measured  and 

calculated  values  of  stresses  compare  favourably  for  the 
two  tests  GC8  and  GC9  which  showed  poor  agreement  using 
the  former  method. 


Using  this  stress  dependent  modulus  for  the  subgrade,  a 
revised  nomograph  was  prepared  and  the  layer  moduli  for  the 
tests  were  re- interpreted.  The  results  of  this,  given  in 
Table  8.3,  show  an  increase  in  gravel  moduli  of  10  -  12%  and 
a  decrease  in  moduli  of  the  subgrade. 

While  we  emphasise  that  the  gravel  testing  program  is  only 
in  the  initial  stage,  we  feel  that  the  results  to  date  are 
useful.  The  series  of  tests  would  benefit  from  additional 
instrumentation . 
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CHAPTER  9 

« 


OPTIMISATION  IN  PAVEMENT  DESIGN 


9.1  INTRODUCTION 


Maximum  economy  in  road  design  can  be  achieved  only  when 
information  on  the  cost  related  performance  of  different 
road  structures  is  available. 

A  systems  analysis  approach  must  be  taken  towards  the 
problem  of  road  design.  This  involves  a  consideration 
not  only  of  the  set  of  possible  initial  designs,  but  also 
the  future  maintenance  requirements  as  well  as  the  possible 
stage  construction  alternatives.  Because  of  the  complexity 
of  this  analysis  and  the  extremely  large  range  of  alternatives 
to  be  considered,  some  form  of  computer  analysis  is  required. 

Two  computer  subroutines,  namely  SELPAV  and  OPTPAV  have  been 
prepared  towards  this  end  in  Trinity  College.  The  former 
is  incorporated  into  the  proposed  design  system  based  on 
the  finite  element  program  DEFPAV ,  (and  is  a  modified  version 
of  that  described  in  Kirwan  and  Maher  (31)).  The  latter  is 
intended  as  an  independent  package  which  considers  the  cost 
implications  of  a  wide  range  of  design  alternatives. 


9.2  ROUTINE  SELPAV 


A  simplified  flow  chart  for  the  routine  is  shown  in  Fig.  9.1. 
It  is  based  on  the  estimation  of  five  critical  parameters. 
These  are  total  construction  thickness,  subgrade  strain, 
tensile  strain  at  the  base  of  the  asphalt  layers,  the 
vertical  stress  on  the  subgrade  and  the  surface  resilient 
deflection . 
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The  search  routine  is  univariant.  Only  the  thickness  of  a 
single  layer  is  varied  at  each  stage.  The  user  specifies  the 
range  of  layer  thicknesses  to  be  considered  as  well  as  the  layer 
thickness  increment.  The  program  is  written  for  the  analysis 
of  a  four-layer  system.  A  three-layer  pavement  can  be  con¬ 
sidered  by  specifying  similar  properties  for  the  two  top  layers. 

The  initial  trial  section  is  the  maximun  possible,  i.e.  that 
comprising  the  maximum  thickness  of  each  layer.  This  con¬ 
figuration  is  then  analysed  to  determine  whether  it  lies  within 
the  allowable  constraints  of  the  five  specified  design  criteria. 

If  the  cross-section  is  found  to  provide  a  workable  solution, 
it  is  costed  and  stored.  The  next  section  is  then  tested. 

The  possible  configurations  are  obtained  by  first  decrementing 
the  thickness  of  the  surface  layer,  until  its  minimum  thick¬ 
ness  is  reached.  Its  thickness  is  then  reset  to  its  maximum 
value  and  the  second  layer  is  reduced  in  thickness  by  the 
specified  decrement.  The  process  of  decrementing  the  surface 
layer  is  then  repeated,  with  each  trial  section  being  tested 
in  turn.  3y  continuing  this  procedure  every  possible  combination 
of  layer  thickness  can  be  tested.  For  a  four-layer  pavement 
with  each  construction  layer  divided  into  ten  increments,  1000 
trial  sections  are  generated. 

To  speed  the  search  process,  if  a  particular  section  fails,  all 
thinner  sections  on  that  cycle  are  not  tested.  To  limit  storage, 
only  the  first  five  hundred  feasible  solutions  are  stored  and 
sorted.  Thus,  if  a  greater  number  of  feasible  solutions  are 
found  (as  specified  in  the  output) ,  the  program  should  be  re¬ 
run  with  larger  layer  increments. 

Estimation  of  Design  Criteria 
Test  1:  Minimum  Pavement  Thickness 

In  this  country,  a  minimum  thickness  of  450mm  of  frost- 
resistant  material  is  recommended  over  frost-susceptible 
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soils.  This  can  be  incorporated,  omitted  or  altered. 


Test  2:  Vertical  Compressive  Strain  in  the  Subgrade 

The  method  estimating  stresses  and  strains  within  the  layered 
structure  is  based  on  the  equivalent  layer  theory  originally 
devised  by  Odemark  (32).  Analysis  of  measured  subgrade  deflections, 
carried  out  by  Ontario  Ministry  of  Transportation  using  a 
deflection  equation  based  on  the  theory,  showed  it  to  be  very 
accurate.  It  has  been  incorporated  into  overlay  design  procedures. 

The  multi-layered  elastic  systerr  is  converted  to  a  homogeneous 
half-space  with  similar  properties  to  the  subgrade.  Thus  the  re¬ 
quired  stresses  and  strains  occurring  in  the  original  layered 
system,  now  occur  in  the  half-space  at  a  depth  related  by  the 
layer  equivalent  thicknesses. 


The  principle  behind  the  calculation  of  the  equivalent  thickness 
he,  of  a  layer  with  thickness  h^,  modulus  E1  and  Poisson's  Ratio 
vt,  in  relation  to  the  material  with  modulus  E 2  and  Poisson's 
Ratio  v 2 ,  is  that  the  equivalent  layer  must  have  the  same  overall 
stiffness  as  the  original  layer,  so  as  to  give  the  same  pressure 
distribution  beneath  it.  The  equations  are  based  on  the  relation¬ 
ship  for  the  modulus  of  rigidity  of  flat  plates  in  flexure. 

E  h3 

D  =  - 

12  (1-v2) 

The  basic  relationship  for  a  two-layer  system  takes  the  form: 
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where  n  is  a.  correction  factor  and  varies  depending  on  the 
number  of  layers  in  the  system.  If  Poisson's  Ratio  is  similar 
for  both  materials,  the  expression  reduces  to: 
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Based  on  the  equivalent  layer  theory,  Ullidtz  (33)  has  formulated 
simple  equations  for  the  estimation  of  stresses,  strains  and 
deflections  in  multi-layered  systems  which  have  given  results 
comparable  to  those  obtained  using  the  computer  program  CHEVRON. 
The  equation  for  subgrade  strain  is  given  by: 


z 


P  (1+  V) 
E 


( 1-2  V)  x 


-  1 


1  +  f-S. 


Where  h  is  the  equivalent  depth  to  the  subgrade  surface,  P  is 
the  applied  stress,  a  is  the  load  radius,  E  is  the  subgrade 
modulus  and  v  the  subgrade  Poisson's  Ratio. 


The  subgrade  strains  for  two,  three  and  four-layer  systems  were 
estimated  using  this  simplified  approach.  The  results  were 
compared  with  those  obtained  using  DEFPAV .  A  range  of  values 
of  modulus,  Poisson's  Ratio,  layer  thicknesses  and  applied 
stress  were  used.  Fig.  9.2  illustrates  the  correlation.  It 
is  considered  satisfactory. 


Test  3:  Tensile  Strain  at  the  Base  of  the  Bound  Layers 

The  equivalent  layer  theory  is  again  used  here.  In  this  case, 
the  subgrade  is  not  taken  into  account  and  so  the  bound  layers 
are  converted  to  an  equivalent  thickness  of  the  sub-base 
material. 

The  vertical  strain  at  the  surface  of  the  sub-base  layer  is 
then  estimated  using  the  equation  for  Test  2.  The  vertical 
stress  at  the  interface  is  found  using  Boussinesq's  equation 
in  terms  of  the  equivalent  height. 
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The  horizontial  tensile  strain  is  estimated  using  the  equation 
(given  by  Ullidtz  (32)  ). 


2  v  ,  (  c  -  sb  e  ) 

sb  '  z  z  ' 


vsb  °z 


£ 

where  sb  is  the  sub-base  elastic  modulus  and  Vsb,  the  sub¬ 
base  Poisson's  Ratio. 


The  accuracy  of  the  method  was  investigated  using  DEFPAV. 

Fig.  4.3  illustrates  the  correlation  for  the  examples  tested. 

Since  the  subgrade  is  not  considered  in  the  equations,  in 
cases  where  there  is  a  large  ratio  between  sub-base  and  sub¬ 
grade  moduli,  larger  errors  occur.  However,  such  modular  ratios 
are  not  common  in  practice  since  they  are  limited  by  the  inability 
of  granular  layers  to  take  tension. 


Test  4:  Vertical  Stress  on  Subgrade 

The  equivalent  layer  theory,  in  conjunction  with  the  Boussinesq 
equation  for  vertical  stress,  is  used  to  determine  the  vertical 
stress  in  the  subgrade.  A  comparison  with  vertical  stresses 
obtained  with  DEFPAV  is  shown  in  Fig.  9.4. 

The  results  can  be  seen  to  be  sufficiently  accurate  for  use  in 
the  selection  routine. 


Test  5:  Surface  Deflections 

A  simplified  deflection  equation,  in  conjunction  with  the 
equivalent  layer  theory,  was  found  to  be  insufficiently 
accurate  for  estimating  surface  deflections  in  a  four-layer 
system.  Therefore,  another  simplified  approach  was  required. 
A  method  based  on  obtaining  an  equivalent  or  mean  modulus  for 
the  multi-layered  system  was  used. 
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The  approach  was  suggested  by  Proksch  (34)  who  used  an 
equation  of  the  form: 


x  hj  +  Ej  x  hj  + 


mean 


hl  +  h2  + 


for  converting  a  multi-layered  system  into  a  two-layer  system 
for  the  purposes  of  estimating  flexural  tensile  stresses  in 
a  design  method.  This  basic  equation  has  been  extended  to 
convert  a  multi-layered  system  into  a  homogeneous  half-space, 
with  an  equivalent  modulus  which  affords  the  same  centre  line 
deflection  as  the  original  system. 


However,  a  problem  arises  when  dealing  with  a  semi-infinite 
subgrade.  It  then  becomes  necessary  to  know  the  depth  of 
pavement  contributing  to  the  the  surface  deflection.  Obviously 
this  will  depend  on  the  thickness  and  elastic  parameters  of  the 
upper  layers  as  well  as  on  the  magnitude  of  the  applied  stress. 
From  preliminary  investigations,  this  influence  depth  was  found 
to  be  approximately  the  depth  at  which  the  vertical  stress 
reached  a  value  of  1/5  of  the  applied  stress.  It  was  also  found 
that,  for  a  system  with  a  constant  applied  stress  and  constant 
material  properties,  the  influence  depth  was  relatively  constant, 
irrespective  of  the  thicknesses  of  the  constituent  layers.  In 
cases  where  the  influence  depth  falls  within  the  upper  layers, 
the  term  referring  to  the  subgrade  is  omitted  and  the  thickness 
of  the  third  layer  is  adjusted.  Thus,  the  equation  for  the 
equivalent  modulus  takes  the  form: 

E1  *  hl  *  E2  X  h2  +  E3  *  h3  +  EsgX  h4 

E  =  - 

eq 

INFLU 

where  E. ,  Z_  and  E,  are  the  successive  layer  moduli,  E  is 
x  di  .j  sg 

the  subgrade  modulus,  and  INFLU  is  the  influence  depth. 

If  h^  +  h2  +  h3  >  INFLU,  h4  =  0  and  h^  is  made  equal  to 

INFLU  -  (hx  +  h2) . 
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The  centre-line  deflection  is  found  using  Boussinesq's 
equation . 


6c 


2  x  P  x  a 


eq 


(1- 


The  variation  of  Poisson's  Ratio  between  the  layers  is  ignored. 


Fig. 9. 5  illustrates  the  accuracy  of  the  results  when  compared 
with  those  obtained  with  DEFPAV. 


The  program  estimates  the  appropriate  influence  depth  given 
the  centre-line  deflection  from  one  trial  section.  Therefore 
a  trial  section  must  first  be  analysed  using  DEFPAV.  That 
consisting  of  the  minimum  layer  thicknesses  is  chosen  for 
convenience.  Since  SELPAV  is  used  in  conjunction  with  DEFPAV 
and  although  surface  deflections  may  not  be  required,  this 
procedure  is  not  considered  to  be  a  disadvantage. 

It  should  be  mentioned  that  the  deflection  estimated  by  the 
routine  is  the  maximum  surface  deflection  under  the  specified 
load.  It  is  not  therefore  equivalent  to  the  Benkelman  Beam 
deflection  which  is  measured  between  two  truck  tyres.  However, 
it  does  permit  a  comparative  assessment  of  the  relative  structural 
capacity  of  the  trial  sections. 


Example  of  SELPAV  Output 

The  input  requirements  are  illustrated  on  the  user  guide 
reproduced  in  Fig. 9. 6.  It  can  be  seen  that  the  user  has  complete 
control  over  the  system  to  be  investigated  and  the  values  of  the 
criteria  used  for  the  final  selection. 

Since  the  analysis  of  a  typical  system  only  takes  a  matter  of 
seconds  of  computer  time,  the  routine  can  be  re-run  many  times 
to  investigate  the  economic  viability  of  alternative  materials 
and  to  determine  which  is  the  critical  parameter  governing  design. 

An  example  of  the  SELPAV  output  for  a  typical  design  problem  is 
shown  in  Figs.  9.7  to  9.9.  The  subgrade  strain  and  asphalt 
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strain  limits  are  those  recommended  by  Shell  (35)  for  a  road 
designed  to  take  a  million  standard  axles. 


ROUTINE  OPTPAV 


Introduction 

The  main  concern  of  a  road-making  authority  is  to  provide  a 
serviceable  pavement  network  at  the  minimum  longrun  cost.  This 
cost  consists  of  the  initial  construction  cost,  maintenance 
costs  and  the  expenditure  on  future  overlays.  For  any  given 
project,  a  large  number  of  alternatives  present  themselves, 
not  only  from  'the  point  of  view  of  choice  of  construction 
materials  and  layer  thicknesses,  but  also  the  initial  design 
life.  In  deciding  on  this,  the  following  points  should  be 
considered. 

1)  Expenditure  on  excessively  long  design  life  pre-empts 
the  use  of  money  for  more  urgent  present-day  projects. 

2)  Pavements  may  become  obsolete  because  of  required  changes 
in  alignment  etc. 

3)  Pavements  rarely  deteriorate  uniformly  along  their  length. 
Thus  most  efficient  use  can  be  made  of  the  available 
budget  by  allowing  for  regular  rehabilitation  on  the  basis 
of  actual  performance. 

4)  The  use  of  shorter  design  lives  allows  greater  flexibility 
of  the  road  management  system,  since  terminal  service¬ 
ability  standards  can  be  revised  to  reflect  variations  in 
available  resources  or  the  relative  importance  of  individual 
roads . 

Hence,  planned  stage  construction  should  offer  the  most 
economical  solution  to  providing  a  serviceable  road  network. 

An  investigation  at  the  T.R.R.L.  concluded  that  two-stage  con¬ 
struction  offered  considerable  advantages  over  initial  design 
periods  of  greater  than  twenty  years. 
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An  optimization  procedure  has  been  devised  which  incorporates 
the  philosophy  of  stage  construction  with  maintenance. 


Basic  Assumptions 

Certain  assumptions  and  simplifications  must  be  made.  In 
order  to  estimate  overlay  requirements,  the  assumption  is  made 
that  the  elastic  properties  of  the  paving  materials  remain 
constant  through  the  pavement  life  and  that  the  effect  of 
cumulative  damage  can  be  taken  into  account  simply  by  its 
effect  on  the  fatigue  life  of  the  bitum  n-bound  layers. 

From  the  literature,  this  assumption  seems  acceptable,  provided 
that  overlaying  takes  place  before  failure  of  any  of  the  pave¬ 
ment  components  has  occurred. 

Burt  (36),  in  a  review  of  progress  in  pavement  design,  concluded 
that  the  properties  of  a  pavement  as  a  whole  do  not  change  much 
during  its  life.  From  an  analysis  of  the  Alconbury  Hill 
experiment  (37) ,  it  was  found  that  the  curve  of  permanent 
deformation  against  cumulative  standard  axles,  after  an  initial 
sharp  fall,  remained  approximately  linear.  While  the  rate  of 
deformation  varied  between  summer  and  winter  and  between  a  good 
and  bad  year  for  weather,  the  average  rate  over  several  years 
was  approximately  constant,  indicating  that  the  material 
properties  were  remaining  constant. 

Measured  transient  vertical  deflections  on  the  same  test  road 
(38)  showed  little  change  with  time  and  traffic  until  after 
critical  conditions  developed.  This  effect  is  also  evident 
from  an  examination  of  overlay  design  curves  based  on  deflect¬ 
ion  (39). 

The  road  management  strategy,  based  on  the  assumption  of 
constant  material  properties,  requires  that  strengthening  is 
undertaken  before  the  onset  of  marked  visible  signs  of 
deterioration.  This  condition  is  the  prerequisite  of  any 
economic  design  approach.  The  design  philosophy  embodied  in 
the  program  OPTPAV  is  an  attempt  to  fulfil  these  requirements. 
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A  simplified  flow  chart  for  the  procedure  is  given  in  Fig. 9. 10. 
The  input  requirements  for  the  program  are  summarised  in  the 
user  guide  shown  in  Fig. 9.11. 


Design  of  Surface-Dressed  Granular  Road 

The  program  initially  considers  a  surface-dressed  granular 
road  to  provide  a  pavement  for  the  minimum  life  specified. 

The  cross-section  is  analysed  using  the  equivalent  layer 
theory,  described  previously  and  the  subgrade  strain 
estimated . 

The  two  relationships  which  relate  allowable  subgrade  strain 
to  pavement  life  in  terms  of  standard  axles  used  in  the 
program  are  those  of  Shell  (35)  and  Brown  et  al  (9).  First, 
using  growth  rate  predictions,  the  number  of  standard  axles 
is  estimated  for  the  trial  gravel  section.  The  subgrade  strain 
is  then  checked  for  the  initial  gravel  section;  if  the  strain 
is  too  large,  the  thickness  is  incremented  until  the  subgrade 
strain  criterion  is  satisfied.  The  subgrade  stress  may  also 
be  specified  as  a  design  criterion. 

When  an  adequate  granular  road  has  been  found,  the  cost  of 
construction,  including  that  of  surface  dressing,  is  estimated. 

It  is  also  necessary  to  estimate  the  maintenance  costs  incurred 
prior  to  the  first  overlay.  This  proves  to  be  a  complex 
commodity  to  estimate  as  shown  by  studies  carried  out  by  the 
Road  Research  Laboratory  (40).  It  includes  such  items  as 
surface  treatment  to  restore  skid  resistance,  renewal  of  road 
markings,  resurfacing  the  hard  shoulder,  raising  kerbs  and 
gulleys  and  local  patching.  Traffic  delay  costs,  as  a  result 
of  maintenance  works,  should  also  be  considered. 

To  combine  this  information  in  the  form  of  a  maintenance 
function  would  form  a  major  study  in  itself  particularly  in 
view  of  the  lack  of  documented  material  on  maintenance  costs 
on  Irish  roads.  As  an  interim  solution,  a  simple  linear 
function  increasing  with  time  is  assumed.  There  is  some  evidence 
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available  (41)  to  indicate  that  this  is  a  reasonable  approxim¬ 
ation. 

In  order  to  make  valid  economic  comparisons  between  alternative 
strategies,  all  future  expenditure  is  discounted  to  present-day 
cost  using  the  equation: 

DRATE  YEAR 

present-day  cost  =  future  cost  x  (1  +  ) 

where  DRATE  is  the  discount  rate  and  YEAR  is  the  number  of 
years  from  present. 

Using  this  relationship,  each  yearly  maintenance  expenditure, 
as  well  as  overlay  costs,  are  discounted  and  summed  to  determine 
total  outlay. 


Design  of  Asphalt  Overlays 

When  the  initial  granular  road  has  been  designed,  the  first 
asphalt  overlay  is  considered.  This  is  designed  on  the  basis 
of  the  allowable  tensile  strain  in  the  asphalt  layer,  the 
allowable  vertical  subgrade  strain  if  desired  and  the  vertical 
stress  on  the  subgrade,  using  the  equivalent  layer  theory  as 
described  previously. 

Then,  using  the  appropriate  fatigue  equation  (four  standard 
equations  at  present  included  in  OPTPAV'  ,  the  allowable  number 
of  standard  axles  are  found.  Miner's  Law  for  cumulative  damage 
is  assumed  to  apply  for  fatigue  failure.  For  the  purpose  of 
the  design  method  in  OPTPAV,  the  pavement  must  undergo  successive 
overlays  (up  to  a  maximum  of  four  in  the  present  version)  with¬ 
out  reaching  a  failure  condition.  Thus  the  concept  of  a 
percentage  of  the  remaining  fatigue  life  to  be  used  for  each 
overlay,  is  introduced. 

If  this  percentage  is  denoted  by  F  (in  the  form  of  a  fraction 
i.e.  F  < 1) ,  then  if  the  first  overlay  has  a  fatigue  life  of 
Nf^  repetitions,  it  is  overlaid  after  reaching  F  x  Nf^ 
repetitions.  The  overlaid  pavement  then  has  a  remaining  life 

1  Of  (1  -  F)  x  Nf 2 . 
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This  in  turn  is  overlaid  when  F  x  (1  -  F)  x  Nf^  repetitions 
have  been  reached  and  continued  for  the  remaining  overlays. 

The  value  chosen  for  F  is  decided  by  the  user.  A  range  of 
values  can  be  tried  to  determine  the  one  yielding  the  most 
economic  result. 

It  is  possible  that  this  theory  of  cumulative  fatigue  damage 
may  produce  too  conservative  a  design  for  future  overlays. 

In  some  existing  design  procedures,  providing  the  facility  of 
stage  construction,  this  damage  theory  is  not  taken  into 
account.  Since,  in  this  country,  fatigue  is  not  the  predominant 
failure  mode  and  also,  since  fatigue  laws,  by  being  based  on 
back  analysis  of  successful  pavement  designs,  have  a  built-in 
factor  of  safety,  ignoring  the  cumulative  fatigue  damage  might 
be  a  legitimate  assumption.  Again,  the  option  of  doing  this 
is  left  to  the  user. 

When  the  first  overlay  has  been  found,  it  is  costed  on  the 
basis  of  construction  and  maintenance  costs  to  the  end  of  the 
design  period.  A  maintenance  function  in  the  same  form  as 
that  for  the  granular  road  is  used.  All  costs  are  reduced  to 
present-day  costs.  The  second,  third  and  fourth  overlay  (if 
required)  are  then  designed  on  a  similar  basis,  until  the 
total  design  life  specified  in  the  input  has  been  reached.  Its 
construction  and  cost  are  then  stored. 

Maintaining  the  initial  granular  road  constant,  the  first 
overlay  thickness  is  increased  from  that  found  in  the  previous 
solution  and  again  the  successive  overlay  requirements 
estimated  and  costed.  The  process  is  continued  until  the 
total  design  life  can  be  reached  by  the  application  of  one 
asphalt  overlay.  Then  the  initial  granular  road  life  is  incre¬ 
mented  and  the  range  of  possible  overlays  determined  as  before. 

When  either  the  initial  granular  road  life  or  its  thickness 
reach  the  maximum  values,  its  thickness  is  reduced  to  the 
minimum.  This  then  forms  the  sub-base  layer  for  the  initial 
three-layer  design  in  which  the  surface  layer  consists  of  the 
minimum  asphalt  thickness.  The  range  of  possible  overlay 


strategies  are  investigated  as  before.  Eventually,  the  sub¬ 
base  thickness  is  incremented  and  the  process  continued. 

In  this  way  every  possible  combination  of  initial  road 
construction  and  life  and  all  possible  overlay  alternatives, 
within  the  limits  specified  by  the  user,  are  found  and  costed. 

On  completion  of  the  search  routine,  the  feasible  solutions 
are  sorted  on  order  of  lowest  overall  cost.  The  first  twenty 
solutions  are  printed  out. 

Figures  9.12  to  9.14  illustrate  the  input  printout  for  a 
typical  design  problem  along  with  a  specimen  output. 
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FLOW  DIAGRAM  FOR  ROUTINE  SELPAV 


Fig.  9.1 


ASPHALT  TENSILE  STRAIN  (DEFPAV)  -  MICROSTRAIN 


L5 


SURFACE  DEFLECTION  (DEFPAV)  -  mm  Q  SUBGRADE  STRESS  (DEFPAV) 


TWO  LAYER 
THREE  LAYER 
FOUR  LAYER 


SUBGRADE  STRESS  (EQUIV.  LAYER)  -  KN/m 

OF  EQUIVALENT  LAYER  THEORY  TO  ESTIMATE  SUBGRADE  STRESS 


0.25 


u  0.25 

SURFACE  DEFLECTION  (EQUIV.  MOD.)  -mm 

ABILITY  OF  EQUIVALENT  MODULUS  THEORY  TO  ESTIMATE  SURFAC 


DEFLECTIONS  Fig.  9 .5 


FORM  OF  INPUT  FOR  ROUTINE  SELPAV 


LOADING  CHARACTERISTICS 


APPLIED  LOAD:-  0.49  8  MN/SQ .  m 
LOADED  RADIUS:-  150.0  mm 

DESIGN  CRITERIA: 


MINIMUM  PAVEMENT  THICKNESS:-  350  mm 
MAXIMUM  TENSILE  STRAIN  IN  BASE  OF  ASPHALT:-  .000150 
MAXIMUM  VERTICAL  STRAIN  IN  SUBGRADE:-  .000650 
MAXIMUM  SUBGRADE  STRESS:-  .0080  MN/SQ .m 
MAXIMUM  RESILIENT  DEFORMATION:-  0.40  mm 


LOADING  AND  DESIGN  INPUT  FOR  ROUTINE  SFLPAV 


LAYER  MATERIAL  CHARACTERISTICS: 


LAYER  1: 

MINIMUM  THICKNESS:-  25  MM 
MAXIMUM  THICKNESS:-  150  MM 
YOUNG'S  MODULUS:-  4000.0  MN/SQ.M 
POISSON'S  RATIO:-  .40 
COST  PER  KM  PER  MM  THICK:-  £  470.00 
LAYER  THICKNESS  INCREMENT:-  25  MM 

LAYER  2: 

MINIMUM  THICKNESS:-  75  MM 
MAXIMUM  THICKNESS:-  300  MM 


YOUNG'S  MODULUS:- 
POISSON'S  RATIO: - 


3000.0  MN/SQ.M 


COST  PER  KM  PER  MM  THICK:-  E  382.80 
LAYER  THICKNESS  INCREMENT:-  40  MM 

LAYER  3: 

MINIMUM  THICKNESS:-  100  MM 
MAXIMUM  THICKNESS:-  700  MM 


YOUNG'S  MODULUS :- 
POISSON'S  RATIO: - 


700.0  MN/SQ.M 


COST  PER  KM  PER  MM  THICK:-  E  30.00 
LAYER  THICKNESS  INCREMENT:-  100  MM 


SUBGRADE : 


AVERAGE  YOUNG'S  MODULUS :- 
POISSON'S  RATIO:-  .42 


30.0  MN/SQ.M 


INFLUENCE  DEPTH  FOR  DEFLECTION  ANALYSIS:-  1907  MM 


INPUT  MATERIAL  CHARACTERISTICS  FOR  ROUTINE  SELPAV 


Fig.  9.1 
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Fig.  9.9 


INPUT 

MATERIAL  PROPERTIES 
TRAFFIC  DATA 
TOTAL  DESIGN  LIFE 
MATERIALS  AND  MAINTENANCE  COSTS 
LAYER  THICKNESS  LIMITS 
MINIMUM  OVERLAY  LIFE 


CHOOSE  INITIAL  GRAVEL  ROAD  LIFE 


DETERMINE  REQUIRED  T 


insmmmi 


ESTIMATE  CAPITAL  & 

I  MAINTENANCE  COSTS  I 

ADD  MIN.  ASPHALT  OVERLAY 


DETERMINE  LIFE  OF  PAVEMENT} 


INCREMENT 
GRA  VEL 
ROAD  LIFE 


INCREASE 

FIRST 

OVERLAY 

THICKNESS 


^  EXTENSION 
> 

s MIN  OVERLAY , 
LIFE. s' 


^ - -  IFF 

TOTAL  DESIGN 
LIFE  ^ 


DETERMINE  FURTHER 
OVERLAY  THICKNESSES 

_ i 

COST 

-  STORE  RESULT 

I 

SORT  IN  ORDER 
OF  LOWEST  OVERALL 
COST 


OUTPUT 


FLOW  DIAGRAM  FOR  OPTRA  V 


INCREMENT 

ASPHALT 

■THICKNESS 


CPA VEL 
THICKNESS  > 
h ALLOWABLE ^ 


SET  GRAVEL 
THICKNESS 
TO  MIN 

y 

SET  ASPHALT 
THICKNESS 
TO  MIN. 


NO  GRAVEL 
ROAD  LIFE 
CONSIDERED 


i 


r.  wo 
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SM 

SP 

Stress 

JSC 

SM:  subgrade  modulus 

SP:  subgrade  Poisson's  Ratio 

Stress:  allowable  subgrade  stress 

JSC:  subgrade  strain  code  1.  Shell  2.  Nottingham 

Default:  1. 


GM 

GP 

GC 

GD 

Q 

R 

Gmin 

Gmax 

Ginc 

GRL 

LMG 

LIFINC 

GM: 

granular  modulus 

GP: 

granular  Poisson's  Ratio 

GC: 

cost  per  km  length  per  mm  thick  as  Layed 

GD: 

cost  of  surface  dressing  per  km  length 

Q.R: 

maintenance  function,  Q  +  R  (AGE) 

G  min: 

minimum  granular  thickness 

G  max: 

maximum  granular  thickness 

G  inc: 

granular  layer  increment 

GRL: 

minimum  granular  road  life 

LMG: 

maximum  granular  road  Life 

LIFINC: 

granular  road  life  increment 

Default:  150  mm 
Default:  3000  mm 
Default:  150  mm 
Default:  1  year 
Default:  15  years 
Default:  1  year 


AM 

AP 

AC 

W 

X 

ASmin 

ASmax 

A  Sine 

MINLIF 

KFC 

AM: 

AP: 

AC: 

W,X: 

AS  min: 
AS  max: 
AS  inc: 
MINLIF: 
KFC: 


asphalt  modulus 

asphaLt  Poisson 's  Ratio 

cost  per  km  length  per  mm  thick  as  layed 

maintenance  function,  W  +  X  (AGE) 

minimum  overlay  thickness 

maximum  total  asphalt  thickness 

asphalt  overlay  increment 

minimum  life  for  asphaLt  overlay 

fatigue  law  code  1.  Shell 


Default:  50  mm 
Default:  1000  mm 
Default:  25  mm 
Default:  5  years 
2.  HR  A 


3.  DBM  (100  pen)  4.  DBM  (200  pen) 


NSD 

Grate 

TDL 

FP 

DF 

Drate 

_ 

NSD: 

GRATE: 

TDL: 

FP: 

DF: 

DRATE: 


number  of  standard  axles  per  day 
percentage  traffic  growth  rate 
total  design  life 

fatigue  life  percentage  Default:  no  cumulative  fatigue 

minimum  thickness  to  prevent  frost  damage  damage  considered, 

discount  rate 


INPUT  USER  GUIDE  FOR  PROGRAM  OPTPAV 


Fig.  9Jj 
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SUBGRADE  CHARACTERISTICS 

RESILIENT  MODULUS:-  70.0  MN/SQ.M 
POISSON'S  RATIO:-  .42 

MAXIMUM  ALLOWABLE  SUBGRADE  STRESS:-  .010  MN/SQ.M 
NOTTINGHAM  CRITERION  FOR  ALLOWABLE  SUBGRADE  STRAIN 


GRANULAR  LAYER  CHARACTERISTICS 

RESILIENT  MODULUS:-  400.0  MN/SQ.M 
POISSON'S  RATIO:-  .35 

COST  PER  KM  LENGTH  PER  MM  THICK  AS  LAYED:-  £  30.00 
COST  OF  SURFACE  DRESSING  PER  KM  LENGTH:-  £  600.00 
MAINTENANCE  FUNCTION:-  £  75.0+40.00  (AGE) 

GRANULAR  LAYER  VARIES  FROM  120.0  MM  TO  1200.0  MM  IN 
INCREMENTS  OF  50.0  MM 

SURFACE  DRESSED  GRANULAR  ROAD  LIFE  VARIES  FROM  3  YEARS 
TO  8  YEARS  IN  INCREMENTS  OF  1  YEAR ( S ) 


SUBGRADE  AND  GRANULAR  INPUT  CHARACTERISTICS  FOR  PROGRAM  OPTPAV 


1 


ASPHALT  LAYER  CHARACTERISTICS 

RESILIENT  MODULUS:-  3500.0  MN/SQ.M 
POISSON'S  RATIO:-  .40 

COST  PER  KM  LENGTH  PER  MM  THICK  AS  LAYED:-  £  800.00 
MAINTENANCE  FUNCTION:-  £  60.0+30.00  (AGE) 

MINIMUM  ASPHALT  OVERLAY  THICKNESS:-  50.0  MM 
MAXIMUM  ASPHALT  THICKNESS:-  300.0  MM 
ASPHALT  OVERLAY  INCREMENT:-  25.0  MM 
MINIMUM  LIFE  FOR  ASPHALT  OVERLAY:-  4  YEARS 
FATIGUE  LIFE  FOR  HOT  ROLLED  ASPHALT  (NOTTINGHAM) 

DESIGN  INFORMATION 


INITIAL  NUMBER  OF  STANDARD  AXLES  PER  DAY:-  600 
TRAFFIC  GROWTH  RATE:-  2.  % 

TOTAL  DESIGN  LIFE:-  30.  YEARS 
FATIGUE  LIFE  PERCENTAGE:-  0.  % 

MINIMUM  DEPTH  TO  PREVENT  FROST  DAMAGE:-  350.  MM 
DISCOUNT  RATE:-  10.0  % 


***  NOTES  *** 

STRESS  AND  MODULI  IN  MN/SQ.M 
MAINTENANCE  ONLY  ESTIMATED  TO  END  OF  DESIGN  PERIOD 
ALL  COSTS  DISCOUNTED  TO  PRESENT  DAY  COSTS 
COSTS  ARE  ESTIMATED  PER  KM  OF  CARRIAGEWAY  10  M  WIDE 


ASPHALT  CHARACTERISTICS  AND  DESIGN  INFORMATION  FOR  INPUT  TO 

PROGRAM  OPTPAV 
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THICKNESS (MM)  1120.  1120.  1020.  1070.  970.  870.  970.  870.  820 

LIFE (YEARS)  6.  3.  0.  0.  0.  0.  0.  0.  0 

CAPITAL  COST ( £ )  34200.  34200.  30600.  32100.  29100.  26100.  29700.  26700.  25200 

MAINTENANCE  COST(£)  885.  378.  0.  0.  0.  0.  0.  0.  0 
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SAMPLE  OUTPUT  FROM  PROGRAM  OPTPAV 
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APPENDIX  1 


ADDITIONAL  INPUT  PARAMETERS  FOR 
ANISOTROPIC  VERSION  OF  DEFPAV 


Vertical  Resilient  Modulus 

The  ratio  of  the  Horizontal  Resilient  Modulus 
to  the  Vertical  Resilient  Modulus, 

The  ratio  of  horizontal  strains  resulting  from 
vertical  strains 

The  ratio  of  horizontal  strains  resulting  from 
horizontal  strains  in  the  orthogonal  direction 

The  ratio  of  Shear  Modulus  in  the  vertical  plane 
to  Vertical  Resilient  Modulus,  G_ 
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